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Abstract 

We present the calculation of the atmospheric neutrino fluxes with an interaction model named 



JAM, which is used in PHITS (Particle and Heavy-Ion Transport code System) 



l|]. The JAM 



^ . interaction model agrees with the HARP experiment [3] a little better than DPMJET-III 

After some modifications, it reproduces the muon flux below 1 GeV/c at balloon altitudes better 
than the modified-DPMJET-III which we used for the calculation of atmospheric neutrino flux in 
5(1 . Some improvements in the calculation of atmospheric neutrino fiux are also 



previous works 
reported. 
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I. INTRODUCTION 



culated atmospheric neutrino fluxes using a 3D-scheme, and refined it, in pre- 
6|. In Ref jo] we studied reputable interaction models with the atmospheric 



We then constructed 



We have ca 
vious works 

muon spectrum at balloon altitudes, and selected DPMJET-III 
and used an inclusive interaction code, which has virtually the same secondary production 
spectra as DPMJET-III, but which runs far more rapidly (see appendix El) • However, when 
we compare the high energy muon flux calculated using DPMJET-III with the observed 
fluxes at sea level and mountain altitudes, there are obvious deficiencies in the calculated 
values at high momenta (> 30 GeV/c) 4]. Therefore, we modified the "inclusive DPMJET- 
III" so that the calculation reproduces the observed atmospheric muon spectra with bet- 
ter accuracy at high momenta, and used it in the calculation of the atmospheric neutrino 
fluxes (5|. We called the modified inclusive interaction code the "modified DPMJET-III". 

n 

Below 5 GeV we used the inclusive NUCRIN which is constructed from NUCRIN [7[ with 
some modifications. However, the contributions of the hadronic interaction below 5 GeV to 
the atmospheric neutrino and muon fluxes even at around 1 GeV/c and below are smaller 
than those above 5 GeV. 

With the "modified DPMJET-III", we could reproduce the observed atmospheric muon 
spectra quite well above 1 GeV/c at sea level and mountain altitudes, but the agreement 
below 1 GeV/c is not as good. At balloon altitudes in particular the calculated muon fluxes 
are obviously lower than the observed ones. This seems to be due to the secondary spectra 
of DPMJET-III, which our simple modification procedure (see appendix [A]) can not correct. 
The disagreement was not resolved without violating the agreement at higher momenta. 
Instead of more sophisticated modification method, we searched for an interaction models 
which has a better nature at low energies, and found the interaction model called JAM, 
which is used in PHITS (Particle and Heavy- Ion Transport code System) JAM is useful 
from much lower projectile energies (<< 0.1 GeV) up to ~ 100 GeV. 

Apart from the atmospheric muon flux, the HARP experiment has greatly improved our 
knowledge of the hadronic interaction at low energies 8|]. It is notable that the hadronic 
interactions of protons on thin N2 and O2 targets have been studied in detail by the HARP 
experiment, aiming to improve the calculation of atmospheric neutrino fluxes j2|. They 
found that the commonly used hadronic interaction models including DPMJET-III do not 



reproduce the result of HARP experiment accurately, especially for the vr productions 
We examine JAM also with the HARP data, and find the JAM agrees a little better than 
DPMJET-HI. 

Then, we construct the inclusive interaction code for JAM and studied the atmospheric 
muon flux with it at sea level, at mountain altitudes, and at balloon altitudes. The "inclusive 
JAM" is used for hadronic interactions from 0.2 GeV to 32 GeV, and above these energies, 
the modified DPMJET-HI is used as before. After some modifications for the "inclusive 
JAM", we can reproduce the observed atmospheric muon flux below IGeV/c at sea level 
and at balloon altitudes, but the agreement of the atmospheric muon flux at the mountain 
altitudes becomes worse. As we can carry out a study similar to the previous publication ^ 
between the muons at balloon altitudes and the interaction model relevant to the low energy 
(< 1 GeV) atmospheric neutrinos, we consider the agreement at balloon altitudes to be more 
important, and use it for the calculation of the atmospheric neutrino fluxes. 

The calculation scheme is essentially the same as that used in Ref 61], although we are 
refining it in each work. In this work, we optimized the "virtual detector correction" to 
minimize the statistical errors. Also with the increase of the available computation power, 
we calculate the atmospheric neutrino fluxes in this 3D scheme up to 32 GeV. We find 
that the calculated atmospheric neutrino fluxes show a visible azimuthal dependence below 
32 GeV. 

The change of interaction model and refinement of the calculation scheme notwith- 
standing, the zenith angle dependence of the atmospheric neutrino fluxes and the the 
(z/^ + v^^jivf. + Pg) ratio are very close to the previous calculations. The change of the 
hadronic interaction model mainly results in a change to the absolute value of the atmo- 
spheric neutrino flux below 1 GeV. 



II. INTERACTION MODEL AND CALCULATION OF ATMOSPHERIC NEU- 
TRINO FLUX 

The requirements for the hadronic interaction model used in the calculation of atmo- 
spheric neutrino fluxes are that it reproduce the real hadronic interaction accurately, and 
that the code generates hadronic events very quickly. The first requirement is a general 
requirement for all applications of the hadronic interaction model. 
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We collected the reputable interaction models, and examined them with available acceler- 
ator data. We further examined them with the atmospheric muon flux observed at balloon 



altitudes. As a result, we selected DPMJET-III for our calculations in Ref [6|. Note, 
DPMJET-III is not usable below ~5 GeV in projectile energy, and so we used NUCRIN 
below 5 GeV. 

The second requirement often conflicts with the first one. The computer code of a so- 
phisticated interaction model generally takes a lot of time to generate an interaction event. 
To resolve this conflict, we construct a quick inclusive interaction code, which generates 
virtually the same secondary production spectra as the original interaction model from the 
output of the original code. We explain the construction of the inclusive code in some detail 
in appendix |Al The interaction code used in Ref jo] was such an inclusive code, made from 
the output of original code of DPMJET-III, and was called the "inclusive DPMJET-III" . 
The inclusive code is useful for the calculation of the time averaged quantity, but is not 
usable for the simulation of an event caused by a high energy cosmic ray, such as an air 
shower. 

The use of the inclusive code has the advantage that we can simply modify the secondary 
spectra of the original interaction model without detailed knowledge of the assumptions 
internal to it. We modified the inclusive DPMJET-III in the procedure explained in ap- 
pendix |X1 and we could reproduce the observed atmospheric muon flux above 1 GeV with 
satisfactory accuracy at sea-level and mountain altitude {4]. The modification changes the 
gradient of secondary spectra globally, and is effectively a change of the so called Z-factor 
{Z =< x^''^ > with X = Esecondary/ Erirojectiie) , which descrlbcs the atuiospherlc muon and 



n 

neutrino fluxes at high energies well [91. The power of 1.7 is approximately the spectral index 
of the integral cosmic ray nucleon flux. We called it modified-DPMJET-III, and calculated 
the atmospheric neutrino flux with it jsl. 

With the modification, however, we could not reproduce the observed atmospheric muon 
flux below ~1 GeV/c. The cosmic ray energies relevant to the atmospheric muon at 1 GeV/c 
are around 10 GeV or lower, depending on the observing altitude, and the cosmic ray 
spectrum can no longer be approximated by a single power law there. The structure of the 
secondary spectra, which does not have much influence on the Z-factor, becomes important. 
The disagreement between calculations and observations in the low energy atmospheric muon 
is due to such a structure in the DPMJET-III model at low projectile energies (< 10 GeV). 




FIG. 1. Comparison of interaction models and HARP data. The solid lines stand for JAM, and 
the dashed lines for DPMJET-III. 



One approach to improving the agreement would be developing a more sophisticated and 
complicated modification procedure. However, we have taken an alternative approach and 
looked for an interaction model which has a better low energy behavior than DPMJET-III, 
and have found an interaction model called JAM which is used in the PHITS (Particle and 
Heavy- Ion Transport code System) 

After our previous calculation, the HARP collaboration released a much more detailed 
study of hadronic interactions of protons on thin N2 and O2 targets at 12GeV/c We 
compared the HARP data with the outputs of DPMJET-III and JAM in FiglH In the figure, 
we combine N2 and O2 target data as the data for air target, assuming the composition of 
air to be 78.5% N2 and 21.5% O2. We also assume the inelastic cross section for proton and 
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FIG. 2. Comparison of the muon flux observed at Tsukuba (30 m a.s.l.) and the one calculated 
with the modified- J AM. Dashed lines in right panel are the ratios of former calculations to present 
ones. 



air nuclei interactions to be 300 mb. In the figure, we find that JAM shows slightly better 
agreement with the HARP experiment than DPMJET-III. Also the secondary spectra of 
JAM is higher than that of DPMJET-III at lower energies. 
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FIG. 3. Comparison of the muon flux observed at Mt. Norikura (2770 m a.s.l) and the one calculated 
with the modified- J AM. Dashed lines in right panel are the ratios of former calculations to present 
ones. 



Constructing the inclusive JAM, we calculate the atmospheric muon flux at the three 
different altitudes where the BESS group has carried out experimental studies of the atmo- 



6 



spheric muon fluxes, at Tsukuba (30m a.s.l.) 10|, at Mt. Norikura (2770m a.s.l.) ll|, and 
at Fort Sumner (balloon altitudes) 12|. As the original JAM seems to have an upper limit 
for the projectile energy (< 100 GeV), we use it below 32 GeV. Above that energy we use 
the inclusive DPMJET-111, as the modified DPMJET-111 reproduces the observed muon flux 
data above >1 GeV/c. 
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FIG. 4. Comparison of the muon flux observed at the balloon altitudes (5-2 g/cm^) at Fort Sumner 
and the one calculated with the modified-JAM. Dashed lines in right panel are the ratios of former 
calculations to present ones. 



We show the comparisons of the calculations using the modified JAM with observations 
in Figs. 121 [3l and HI as the ratios of the observed muon fluxes to the calculated ones. In these 
calculations, we used the air density profile measured on the same date as the observations. 
This makes the agreement between the calculation and observation better for >lGeV/c. 
Some modifications are also applied to the JAM model, following the procedure explained 
in appendix |Al 

We note that it was very difficult to modify the inclusive DPMJET-III only to give the 
muon flux as large as the observed one at Tsukuba and at Fort Sumner observed by BESS 
below IGeV/c. Introducing the inclusive JAM with the some modifications, we can boost 
the calculated muon flux below 1 GeV, and get good agreement between the calculations and 
observations at Tsukuba and Fort Sumner. However, boosting the muon flux below 1 GeV 
makes the difference between the calculation and observations at Mt. Norikura larger. We 
consider that it is more important to reduce the difference at balloon altitudes, since the 
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atmospheric muon at the balloon altitude is directly related to the primary cosmic rays and 
the hadronic interactions at the energy relevant to the low energy (< 1 GeV) atmospheric 
neutrinos. We can carry out a study of the hadronic interaction at this energy as we 
did in the the previous publication using the muons observed at ground level at the 
energy relevant to the atmospheric neutrino with energies above 1 GeV. We calculate the 
atmospheric neutrino flux with the modified JAM in this paper. 




FIG. 5. Comparison of pion production spectra between the original JAM (solid line) and the 
original DPMJET-III (dashed line) integrated all over the directions in the left panel, and those 
of the modified-JAM (solid line) and the modified-DPMJET-III (dashed line) in the right panel at 
10 GeV. Here the X denotes the ratio of secondary kinetic energy to projectile kinetic energy. 



In Fig. IHl we show the comparison of the x-spectra (x = Egecondary/ Eprojectue) of the orig- 
inal JAM and the original DPMJET-III, and those of the modified-JAM and the modified- 
DPMJET-III which we actually used to calculate the fluxes in Figs. 121 |3l and HI In both 
panels, we find that JAM shows slightly higher spectra than DPMJET-III in x < 0.2. This 
is the reason for the higher atmospheric muon flux below 1 GeV/c with the modified-JAM. 
And, our modification method cannot modify DPMJET-III to have this high x-spectra for 
X < 0.2 and also keep the high spectra for x > 0.2 required by the observed muon flux above 
1 GeV. 



III. CALCULATION OF ATMOSPHERIC NEUTRINO FLUX 

The scheme for calculating the atmospheric neutrino fluxes is almost the same as for 
the previous calculations j^, 0], and so we only briefly review it here. We use the same 
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of —2.71 above 100 GeV (see also Refs. 



primary flux model based on AMS 13|, [14 [and BESS 10|, ll5| data, with a spectral index 



16 



17|). We used the US-standard '76 model 



of the atmosphere 18|, as the error due to the atmospheric density model is sufficiently 



19| geo- 



small for the calculation of atmospheric neutrino fluxes ^ . We use the IGRF2005 
magnetic model, as the variation due to changes of geo-magnetic fleld would be small, and 
only extrapolation is available to the year 2010 at the time of computation. 

For the 3D calculation of the atmospheric neutrino fluxes, we work in the Cartesian 
coordinate system which has the origin at the center of the Earth, with the Z-axis extending 
to the north pole, and we consider the surface of the Earth is a sphere with a radius of 
Re = 6378.180 km. However, the position on the Earth is rather well represented by the 
spherical polar coordinate system (r, 6, 0) with r = R^, related to the Cartesian coordinate 
system by 

X = Re sin 6 cos (j), 

y = Re sin 6 sin 0, and (1) 
z = Re cos 6 

The local coordinate system at the detector is constructed based on this polar coordinate 
system. The direction of the x-axis is in the increasing direction of 6, the direction of the 
y-axis is in the increasing direction of </>, and the direction of the z-axis is in the increasing 
direction of r. Therefore, the azimuth angle is measured counterclockwise from south in the 
local coordinate system. 

In addition to the surface of the Earth, we assume three more spheres; the injection 
sphere, the simulation sphere, and the escape sphere. In the previous work, we took the 
radius of the injection sphere as Rinj = Re + 100 km, the radius of simulation sphere as 
Rsim = Re + 3000 = 9378.18 km, and the escape sphere as Resc = 10 x Re = 63781.80 km. 
However, we take the radius of the simulation sphere to be Rsim = 10 x Re = 63781.80 km, 
the same as the escape sphere, as a result of the increase in computation power. However, 
this change does not cause any visible difference in the results. 

Cosmic rays are sampled on the injection sphere uniformly toward the inward direction, 
following the given primary cosmic ray spectra. Before they are fed to the simulation code 
for the propagation in air, they are tested to determine whether they pass the rigidity cutoff, 
i.e., the geo-magnetic barrier. For a sampled cosmic ray, the 'history' is examined by solving 
the equation of motion in the negative time direction. When the cosmic ray reaches the 
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escape sphere without touching the injection sphere again in the inverse direction of time, 
the cosmic ray can pass through the magnetic barrier following its trajectory in the normal 
direction of time. 

The propagation of cosmic rays is simulated in the space between the surface of Earth 
and the simulation sphere. When a particle enters the Earth, we discard the particle. Note, 
those particles which enter the Earth are mainly muons, and the number of muons which 
enter the Earth are far smaller than the number of muons which decay in air. Also, when 
muons enter the Earth, they lose energy quickly, and decay at rest or are captured by an 
atom in the rock and produce a neutrino with energies lower than 100 MeV. 

Note, the neutrino detectors are very small compared with the size of the Earth, and 
are considered as infinitesimal points on the surface of the Earth. We need a finite size 
"virtual detector" instead of the infinitesimal neutrino detector in the 3D calculation of the 
atmospheric neutrino flux. However, the averaged neutrino flux in the finite size "virtual 
detector" is a little different from the neutrino flux at the actual neutrino detector, and we 
introduced the ■.vi.tua. detect. eo„.ect,o.- in Ref. fl. TKi. "vutual deteeto. co^eetion- 
is refined in this paper. 

n 

The procedure for the correction is the same as in Ref. [5|. Using the inside of the circle 
with radius of 9d as the virtual detector, the flux $6*^ obtained with the virtual detector and 
the flux $0 at the real target detector are related by 



^ $0 + ^',el , (2) 

The term with $q may be cancelled out using two virtual detectors, with radii of 6i and 621 



as 

32 



°- el -el i-r2 ' ^""^ 

where we assumed 6'i > 6*2. 

In appendix [Bl we study the ratio of the statistical errors A$o/^*^^i the function of r. 
We find that the ratio takes the minimum value v^S) ~ 2.236 at r = 1/^2 ~ 0.707, whereas 
we took r = 0.5 in Ref. jsj]. However, the ratio A^o/^"^! is a slowly varying function of r 
near the minimum, and the ratio at r = 0.5 is not too bad (2.517). 

We take the radius of the virtual detector 1 as 1113.6 km corresponding to the angle 
from the center of the Earth of 10 degrees, as in Ref. jsl, but we take the radius of the 
virtual detector II as 787.4km corresponding to 7.071 degrees in this work. The "virtual 
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detector correction" is applied up to E^, = 10 GeV, as we find a large difference in the fluxes 
determined by the virtual detectors I and II for some directions near the horizon (Fig. [T71 
in appendix [B]) . 

In the previous work, we connected the atmospheric neutrino fluxes calculated in the 3D 
scheme to the one calculated in the ID scheme at 10 GeV. In this work, we connected the 
3D scheme calculation to ID one at 32 GeV, since we find a dipole-like variation of z/g and 
Ufj_ for horizontal directions due to the muon bending in the geo-magnetic field at 10 GeV. 

It is not efficient to calculate the atmospheric neutrino spectrum from 0.1 GeV to 32 GeV 
in one simulation, as the cosmic ray spectra are so steep. Therefore, we calculate the 
atmospheric neutrino flux limiting the lowest neutrino energy to 1 GeV, 3.16 GeV, 10 GeV, 
and 31.6 GeV, also limiting the lowest energy of primary cosmic ray to 1 GeV, 3.16 GeV, 
10 GeV, and 31.6 GeV respectively. The flux calculated in those simulations are combined 
into single atmospheric neutrino spectrum for each direction and each kind of neutrino. 

We show the statistical error in this work for near-vertical directions and for near- 
horizontal directions averaging over all azimuth angles, in Fig. O For the horizontal di- 
rections, the opening angle of the virtual detector decrease with cos 6, where 6 is the zenith 
angle. Therefore, the statistical error for near-horizontal directions is larger than that for 
near-vertical directions. The statistical error depends not only on the calculation scheme, 
but also on the analysis method of the simulation data. The analysis method of the sim- 
ulation data and the estimation of the statistical error in our calculation is explained in 
appendix O in some detail. 

We find the statistical errors are well below 1 % for the average over the azimuth angles. 
For each azimuth bin, however, the statistical errors reach a few percent for some directions. 
The largest statistical errors are found at the energy just below 10 GeV, due to the virtual 
detector correction applied up to this energy. 

IV. ATMOSPHERIC NEUTRINO FLUXES CALCULATED WITH MODIFIED- 
JAM 

We show the calculated atmospheric neutrino fluxes with the modified- JAM in this sec- 
tion. In the left panel of Fig. [TJ the calculated atmospheric neutrino spectra are shown 
averaged over all directions from 0.1 GeV to 32 GeV for Kamioka, together with the ones 
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FIG. 6. The ratio of statistical errors of the atmospheric neutrino flux calculated to the flux value 
in this work. We show the statistical errors for near- vertical directions (0.9 < cos^z < 1) in the 
left panel, and for near-horizontal directions (0.9 < cos^^ < 1) in the right panel, both averaging 
over all azimuthal directions. 



calculated in our previous work, and those of the Bartol group 



group 



HQ, 



and the FLUKA 



22|. Above 32 GeV, the 3D calculation is smoothly connected to the ID calculation 



carried out in the previous work. As the modified- JAM is used below 32 GeV, any difference 
above that is due to the difference of the calculation scheme between 3D and ID. However, 
the difference between present and previous works is very small in the figure above 1 GeV. 
Note, we magnify the difference between present and previous works in the ratio in Fig. [HI 
The difference is less than a few percent above 1 GeV. On the other hand, the atmospheric 
neutrino fluxes calculated with the modified- JAM show an increase from the previous one 
below 1 GeV, as is expected from the increase of atmospheric muon spectra below 1 GeV/c. 

In the right panel of Fig. [TJ we show the ratios of the atmospheric neutrino flux averaged 
over all directions. It is seen that the differences in the flux ratios are very small between 
present and previous works above 1 GeV as the absolute atmospheric neutrino fluxes. It 
is noticeable that the {u^ + z/^) / (z/g + z/g) ratios are very close to each other including the 
ones calculated by different authors. We find there is a visible difference in Ue/i^e between 
the present and previous works. As is seen in the left panel of Fig. [S] in sec [TTl the original 
JAM interaction model has a smaller 7r"*'/7r~ ratio, and is responsible for the smaller z/g/z/e 
ratio below 1 GeV. Note, we do not modify the interaction model more than the muon flux 
data require.. It is difficult to examine the fJ.'^ / fJ.~ ratio at the momenta corresponding to 
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FIG. 7. Comparison of atmospheric neutrino fluxes calculated for Kamioka averaged over all 
directions (left panel), and the flux rati os ( right panel), with other calculations. The dashed lines 



are the calculation by the Bartol group 



2C| 



2l[, dotted lines for the FLUKA group 



22| . and dash 



dot for our previous calculation (HKKM06). 



the Vf^j Vf. below 1 GeV, due to the small statistics in the observation at balloon altitude. 

The muons at sea level or mountain altitude are not useful to examine the atmospheric 
neutrino of these energies, since the muons result from higher energy pions at higher altitude. 

In Fig. [H] we show the atmospheric neutrino fluxes as a function of the zenith angle 
averaging over all the azimuthal angles at 3 neutrino energies; 0.32, 1.0, and 3.2 GeV for 
Kamioka. In Fig [10] we show the comparison of the present and previous works in the ratio 
as the function of zenith angle. There is a difference due to the increase of the neutrino flux 



le calculated zenith angle dependences 
6|- 



itself, but the ratio is almost constant. Actually, t 
are virtually the same as for the calculation in Ref 

It seems that the zenith angle dependence of the 3D calculation smoothly connected 
to that of the ID calculation just above 3.2 GeV for the average over all azimuth angles. 
However, this is not true when we study the variation of atmospheric neutrino flux as a 
function of azimuthal angle. In Fig. [TT] we show the variation of atmospheric neutrino flux 
as the function of the azimuthal angle averaging them over the five zenith angle ranges, 
1 > cos^ > 0.6, 0.6 > cos^ > 0.2, 0.2 > cos^ > -0.2, -0.2 > cos^ > -0.6, and 
—0.6 > cos ^ > — 1, at 1 GeV for Kamioka. It is seen in that the variation of the atmospheric 
neutrino flux has complex structures at 1 GeV due to the rigidity cutoff and muon bending 
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FIG. 8. Atmospheric neutrino fluxes calculated for Kamioka averaged over all directions, as a ratio 



to the HKKM06 calculation. The dashed lines are the calculation by Bartol group 



dotted lines the FLUKA group 



23]. 



and 



in the geo-magnetic field 23|]. Note, the variation of upward going neutrinos (—0.2 > cos^) 
is much more complicated than the variation of downward going neutrinos {cos 6 > 0.2). 
This is because the upward going neutrinos are produced in a far larger area on the Earth 
than the downward going neutrinos, and are affected by large variation of rigidity cutoff and 
geo-magnetic field. The downward going neutrinos are produced just above the detector. 

In Fig. [12] we show the variation of atmospheric neutrino flux as the function of the 
azimuthal angle averaging them over the same five zenith angle bins at 3.2 GeV for Kamioka. 
At this neutrino energy, we find the deficit of neutrino flux due to the rigidity cutoff from 
the East direction (</> ~ 90°) clearly for all kind of neutrinos and all zenith angle bins. 
However, the z/g and z/^ show large dipole structures in the horizontal zenith angle bin 
(0.2 > cos^^ > —0.2) due to the muon bending in the geo-magnetic field. 

In Fig. [13] we show the variation of atmospheric neutrino flux as a function of the az- 
imuthal angle averaging over the same flve zenith angle bins at 10 GeV for Kamioka. At this 
energy, the atmospheric neutrinos seem to be free from the rigidity cutoff, but the z/g and 

still show the dipole structures in the horizontal zenith angle bin (0.2 > cos 6' > —0.2). 
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FIG. 9. The zenith angle dependence of atmospheric neutrino flux averaged over ah azimuthal 
angles calculated for Kamioka. Here 9 is the arrival direction of the neutrino, with cos 9 = 1 for 
vertically downward going neutrinos, and cos 9 = — 1 for vertically upward going neutrinos. 
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FIG. 10. Ratio of the atmospheric neutrino flux averaged over azimuth angles as the function 
of zenith angle to the same quantity calculated in HKKM06. The short dash lines stand for the 



calculation by Bartol group 
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2l[ and dotted lines for the FLUKA group 
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The dipole structure of and seems to extend to a few 10 GeV in the horizontal 
bin. In Fig. [TH we show the variation of atmospheric neutrino flux as the function of the 
azimuthal angle averaging them over the same five zenith angle bins at 32 GeV calculated for 
Kamioka. We still see the dipole azimuth angle variation of z/g at this energy in the horizontal 
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FIG. 11. The azimuthal angle dependence of atmospheric neutrino flux at 1 GeV, averaged over 
zenith angle bins of 1 > cos 6' > 0.6, 0.6 > cos6' > 0.2, 0.2 > cosi9 > -0.2, -0.2 > cos6' > -0.6, 
and —0.6 > cos0 > —1 calculated for Kamioka. 



bin. However, the amplitude of the variation is almost the same as that of the statistical 
error. Therefore, we connect the atmospheric neutrino flux calculated in 3D scheme to the 
ones calculated in ID scheme at 32 GeV. In Fig. [HI we also plotted the flux vales calculated 
in the ID scheme. We can say the flux values calculated in the 3D and ID schemes agree 
within the statistical error, except for the dipole structure of v^. in the horizontal bin. 



V. SUMMARY AND DISCUSSION 

After the previous publication, we have searched for an interaction model which has 
better behavior at low energies than DPMJET-III and NUCRIN, and we found that the 
JAM interaction model, used in PHITS (Particle and Heavy- Ion Transport code System) 
It shows slightly better agreement with the HARP experiment than DPMJET-III. 

Constructing the inclusive interaction model from the JAM interaction model, we have 
used it as the hadronic interaction code in the simulation of cosmic ray propagation through 
the atmosphere below 32 GeV. Above this energy, we continue to use the modified DPMJET- 
III, since we could reproduce the result of atmospheric muon observations by BESS for 
> 1 GeV/c in the previous work. 
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FIG. 12. The azimuthal angle dependence of atmospheric neutrino flux at 3.2 GeV, averaged over 
zenith angle bins of 1 > cos 6* > 0.6, 0.6 > cos6' > 0.2, 0.2 > cosi9 > -0.2, -0.2 > cos6' > -0.6, 
and —0.6 > cosO > —1 calculated for Kamioka. 

Applying some modifications to the inclusive JAM, we were able to increase the atmo- 
spheric muon flux below 1 GeV/c, and achieved better agreement with BESS observations 
at Tsukuba (30m a.s.l.) and at Fort Sumner (balloon altitudes). However, the increase of 
the atmospheric muon flux below IGeV/c made the difference larger between calculation 
and observations at Mt. Norikura (2770m a.s.l.). We consider the agreement at the balloon 
altitudes to be more important, since muons are observed there with almost the same mo- 
mentum as they are produced. We then calculate the atmospheric neutrino fluxes with the 
modified- JAM as the hadronic interaction code below 32 GeV. 

The calculation scheme is essentially the same as in previous calculations, but the radius 
of the simulation sphere has been made the same as that of the escape sphere, to reduce 
the boundary effect. However, there is no visible effect from this change in the calculated 
atmospheric neutrino fluxes. Also we revised the formula for the size correction of the virtual 
detector to minimize the statistical errors, which were reduced by 

Some features of the calculated atmospheric neutrino flux are presented in section IIVI 
The use of the modified- JAM increased the atmospheric neutrino fluxes below 1 GeV corre- 
sponding to the increase of atmospheric muon fluxes below 1 GeV/c. However, the ratios of 
the atmospheric neutrinos, especially (z/^ + P^)/(z/e + Ue), are very close to the previous cal- 
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FIG. 13. The azimuthal angle dependence of atmospheric neutrino flux at 10 GeV, averaged over 
zenith angle bins of 1 > cos 6* > 0.6, 0.6 > cos(9 > 0.2, 0.2 > cos6' > -0.2, -0.2 > cos6i > -0.6 
and —0.6 > cos 6 > —1 calculated for Kamioka. 

culation. Only the ratio i^e/^'e shows a decrease near 0.1 GeV due to the small tt^/tt^ ratio 
of the original JAM. It is difficult to study the vr^/vr" ratio from the observed atmospheric 
muons. Also the zenith angle variations of the atmospheric neutrinos are virtually the same 
as the previous calculations. 

The azimuth angle variations of the atmospheric neutrino fluxes are studied in detail in 
this paper. For horizontal directions (0.2 > cos 6* > —0.2), and z/g still show a dipole-like 
variation even at 10 GeV due to the deflections of muons in the geo-magnetic field. However, 
the magnitude of the variation become almost the same as the statistical error at 32 GeV 
even for horizontal directions. As a result, we connected the atmospheric neutrino fluxes 
calculated in the 3D calculation scheme to those calculated in the ID calculation scheme at 
32 GeV. They agree with each other within the statistical error at this energy. 

The study of the interaction model and atmospheric muons, similar to that in the previous 
publication, suggests that the magnitude of the error of the atmospheric neutrino flux due to 
the hadronic interaction model at energies below 1 GeV is similar to the maximum difference 
between the calculated atmospheric muon fluxes and observed ones below < 3 GeV/c at 
balloon altitudes. Assuming ~10% of the uncertainty for the total interaction cross section 
of cosmic rays and air nuclei, or the interaction mean free path of the cosmic rays, we 
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FIG. 14. The azimuthal angle dependence of atmospheric neutrino flux at 32 GeV, averaged over 
zenith angle bins of 1 > cos 6* > 0.6, 0.6 > cos6' > 0.2, 0.2 > cos6l > -0.2, -0.2 > cos6' > -0.6, 
and —0.6 > cos^ > — 1 calculated for Kamioka. The straight lines in the flgure show the flux 
values calculated in the ID scheme. 



estimate the uncertainty of the atmospheric neutrino flux below 1 GeV is around 15% at 
0.3 GeV, and 20% at 0.1 GeV. 
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Appendix A: Interaction model and calculation of atmospheric neutrino flux 

Here, we explain the construction of the inclusive code from the output of the original 
interaction code. First we ran the original code many times, and recorded the kind, energy, 
and momentum of the secondary particle. We sort and number them in the ascending order 

of — Eggi^Qndary / Eprojectile 

Xi < ■ ■ ■ < Xi-i < Xi < Xi+i < ■ ■ ■ < Xn , (Al) 
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for each kind of secondary particle. Considering this series as a function of integer in [l,n\, 
we can define a function of a discrete variable in (0, 1] as 

X. = X(^) . (A2) 

This X{i/n) converges to a continuous function in [0,1] as the number of trials, or n, 
increases, and the asymptotic continuous function can be approximated by a proper fitting 
function for sufficiently large n. 

Writing the asymptotic function also as X{u), we can generate the x following the x- 
distribution of the original interaction model with uniform random numbers m in [0, 1] as 

x = X{u) . (A3) 

It is straightforward to calculate the energies of the secondary particles from x. 

As the fitting function, we used the B-spline function for DPMJET-III and NUCRIN, but 
we take poly-quadra function for JAM for accuracy, although the difference in the actual 
usage is small. 

We fit the angular distribution by the fitting formula, 

F(cos^)dcos^ = exp(A + a(cos^ + l)"(l-cos=^^)^)dcos^ , (A4) 

at low energies (< 32 GeV). We determine the parameters, A, a, a, and /3 for each x-bin 
with log Ax = 0.1. For higher energies, we fit the angular distribution with a simple formula 

F'(cos 9)d cos 9 = exp{A + a ■ cos 6)d cos 6 , (A5) 

which is a special case of eq. IA4[ 

Repeating the above procedures for all kinds of secondary particle and for all the projectile 
energy grids, the inclusive code is constructed for the interaction model. Note, the projectile 
energy grids are selected to cover the required projectile energy range. For the projectile 
energy in between two energy grids, we use the set of parameters at either grid energy. The 
grid energy to use is determined by the distances to both grid energies and a random number 
in a probabilistic way. 

As the grid projectile energies, we took log(-E'g/l GeV) = 0.5, 1.0 ... 5.5, and 6.0 for 
DPMJET-III, and \og{Eg/l GeV) = -0.7, -0.6, -0.4, .0, and 0.5 for NUCRIN. Fro JAM 
interaction model, we take log{Eg/l GeV) = -1, -0.9, ... 0., 0.1, ...1.6, and 1.7. 
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As noted previously, the inclusive code allows us to readily modify the secondary spectra 
of the original interaction model. When we apply a modification to the inclusive code, we use 
replace the uniform random generator to the weighted random number generator in [0, 1], 
instead of the uniform random number generator. We take a simple form for the weight 
function, 

Wi{u) = 1 + ai ■ u , (A6) 

and request the random numbers are generated in the probability proportional to w{u). The 
Oj are the modification parameter defined at the projectile energy grid for all the kind of 
secondary particles. At each projectile energies, the are not independent but the total 
energy of secondary particle is conserved on average. The multiplicities are not modified. 

The above modification procedure is an effective way to modify so called Z-factor [Z =< 
x^'^ >, with X = Esecondaryl E projectile) which is a good description for atmospheric muons 
and neutrinos at high energies |9(]. 



Appendix B: "Virtual detector correction" and the optimization 



We have shown in Ref. [51] that the average flux ($6^) ^ virtual detector with radius of 
6d is expressed as 

^ $0 + <^',el , (Bi) 

where $o is the atmospheric neutrino flux at the target detector, and $[, is a constant 
determined for the target detector. 

As in Ref. [sl, we may cancel out the term with $q using two averaged fluxes in two 
concentric virtual detectors around the target detector. Assuming two virtual detectors 
with radii 9i and 6'2, we can calculate the flux at the target detector ($o) as 

where, $i is the averaged flux in the virtual detector I with a radius of 6*1, and $2 is the 
averaged flux in the virtual detector II with a radius of O^- Here, we assumed 9i > 62 

The treatment of the virtual detector in Ref. jsj] corresponds to taking r = 0.5 in |3l 
However, this value is not optimal for the statistical error due to fluctuations in the number 
of neutrinos observed in the virtual detector. The fluxes $1 and $2 are calculated as 



^1 . ^ ^: 



$1 = ^ and $ 



2 
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where T is the corresponding observation time of the simulation, and A^i and N2 are the 
observed numbers of neutrinos in the virtual detectors 1 and II respectively. The "true" flux 
$0 is calculated from Ni and as 



$0 



(B3) 



Ttt91 ■ (1 - r2) ■ 

We note that the neutrino observed by the virtual detector II is also observed by virtual 
detector I, therefore, A'^i — N2 and N2 are independent variables rather than A'^i and 
Therefore, in terms of independent variables, the "true" flux $0 is rewritten as. 



*° Tnef ■ (1 -r2^ ^ 



(--r'')N2-r''{N^-N2) 



(B4) 
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FIG. 15. The variation of the function F{r) defined bv lB6[ It has a minimum value ~ 2.236 
at r = 1/v^ ~ 0.707 indicated by "Optimized". The "HKKM06" indicates the r = 0.5 used in 
HKKM06, where the -F-value is a little larger than the optimized one (-F(0.5) = 2.517). 

Assuming Poissonian fluctuations for the numbers of the observed neutrinos, we can 
estimate the statistical error of the "true" flux as 



1 



Tn9l ■ (1 - r2) 



(B5) 



where we used the approximate relations — r'^Ni and Ni — N2 — {1 — r'^)Ni. Therefore 
the ratio of A$o/A$ 1 IS given as 



A$i 1 - r 



1 ''i- -r2)2r2 + r4(l -r2) = F(r) . 



(B6) 



The function F{r) has a minimum at r = 1/^2 ~ 0.707 and ^(1/^2) = v^5) ~ 2.236. We 
take the radius of the virtual detector I as 1113.6 km corresponding to an angle of 10° at 
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the center of the Earth, and the radius of the virtual detector II as 787.4 km corresponding 
to 7.071°. 

The "virtual detector correction" is studied in the [before/after] ratio in our calculation 
for Kamioka in Fig. [T51 for 1 > cos^^ > 0.9 as the near vertical directions and for 0.1 > 
cos ^2 > as the near horizontal directions, where 9^ is the arrival zenith angle. Note, the 
"virtual detector correction" is necessary and effective for the downward going neutrinos . 
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FIG. 16. Ratio of the atmospheric neutrino flux averaged over the virtual detector to 
the flux for the true neutrino detector ((/>(0)) estimated with the procedure explained in the text. 
The atmospheric neutrino fluxes are averaged over all the azimuth angles. Here, 9z stands for the 
zenith angle of the arrival direction, and the left panel is for the near vertical downward directions, 
and the right panel for the near horizontal directions. 

From the figures, we find the "virtual detector correction" is needed below 3 GeV even for 
the average over all azimuth directions. When we divide all the arrival azimuth angles into 
12 bins, we find the largest correction in the azimuth bin of [90°, 120°] for near horizontal 
directions (0.1 > cosOz > 0). Note, we measure the azimuth angle counterclockwise from 
south. We show the [before/after] ratio of the flux for this bin and opposite azimuth bin 
([270°, 300°]) in Fig. [T71 summing over all kinds of neutrino. 

Appendix C: Analysis of simulation data and the statistical error 

To analyze the simulation output, we take a logarithmically equally-spaced grid Ei for 
the neutrino energy as 

Ei = Eo-R\ (z = ..,-1,0,1,....), (CI) 
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FIG. 17. Ratio of the neutrino fluxes averaged over the virtual detector {(j){10)) to the flux for the 
true neutrino detector (0(0)) calculated with the procedure explained in the text for near horizontal 
east (90° < (pa < 120°) directions. 

and count the the numbers of neutrinos falling in each [-Ej+i, Ei] bin for each kind of neutrino 
and each arrival direction window. Note, we take Eq = 0.1 GeV and R = 10^^^", then used 
grids with z > — 1 in this study. As usual, we assume that the integration of the neutrino 
flux in a energy bin is related to the number of neutrino in the bin as 



+ 1 



^{E)dE 



s-n-T 



(C2) 



where Ni is the number of neutrinos in the [Ei+i,Ei\ bin, T is the equivalent time of the 
simulation, Q is the opening angle of the window, and S is the area of the virtual detector. 

We calculate the neutrino flux $j at energy Ei from /j and which are integral fluxes 
calculated from the simulation in the two adjacent energy bins [Ei+i,Ei\ and 
respectively. For these two energy bins, we approximate the neutrino flux as the power 
function of the energy as $(i?) = A ■ E°', and $j is calculated as 

(ln(J,)-ln(J,_i)) 



(C3) 



EMR){i^-h^i) 

This expression is singular at Jj = but it becomes a smooth function of Jj and Ji_i, if 
we define the value at Jj = as 



EMR) 

to which IC3I approaches in the It — )■ limit. 



(C4) 
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Assuming Poissonian fluctuations for the statistical error, the error for integrated flux 
l^lijli is given by Xj^fWi. The statistical error for $j is estimated substituting IC3l into the 
general expression 

and the ratio A$j/$j is given as 



A<l>,- //I 1\M / I 



(C6) 



$i V Vlnx x-\) Ni \\n{l/x) 1/x-lJ iVi_i ' 
where x = li/Ii^i = Ni/Ni_i. This expression is well approximated by the simpler relation, 



In the Ni — > Ni^i limit, both [C6] and ICTl give 1/y/Ni + A^2, and the differences between 
them are less than 5 % for 0.5 < Ni/Ni_i < 2. 



Appendix D: Atmospheric neutrino flux below 32 GeV 

Here we tabulate the atmospheric neutrino flux for solar minimum at sea level at locations 
— Kamioka, Sudbury (North America), and Gran Sasso — as in the former publication, 
averaging them over all the arrival azimuth directions and in the zenith angle bins with 
A cos 6*2 = 0.1, where 6^ is the zenith angle of the arrival direction of the neutrino, in Tables 
I-XX. 

Due to page limitations, we cannot present all the calculated results here. The atmo- 
spheric neutrino fluxes for different solar activity for Kamioka, Sudbury, Soudan2 site. Gran 



Sasso, Frejus, and Homestake will be available at the web site: http : / /www . icrr . u-tokyo . ac . jp/~mhonda 
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TABLE I. Neutrino flux (m-^sec-^sr-^GeV-^) for 1.0 > cos > 0.9 



(GeV) 



Kamioka 



Sudbury 



Gran Sasso 



Norm 



4.25 4.25 2.02 2.08 

3.52 3.52 1.66 1.68 

2.81 2.81 1.31 1.31 

21.50 21.50 9.98 9.81 

15.68 15.68 7.32 7.04 

11.02 11.02 5.13 4.84 

7.46 7.46 3.47 3.21 

4.88 4.88 2.26 2.05 

3.09 3.09 1.42 1.26 

18.93 18.93 8.65 7.54 

11.31 11.31 5.06 4.32 

6.55 6.55 2.87 2.40 

3.70 3.70 1.58 1.29 

20.28 20.28 8.36 6.67 

10.83 10.83 4.25 3.36 

5.68 5.68 2.11 1.64 

29.33 29.33 10.13 7.77 

15.07 15.07 4.76 3.62 

7.65 7.65 2.19 1.66 

38.61 38.61 9.94 7.52 

19.36 19.36 4.43 3.37 

9.70 9.70 1.97 1.50 

48.46 48.46 8.68 6.71 

24.26 24.26 3.84 2.99 

12.07 12.07 1.71 1.33 

60.83 60.83 7.68 6.08 



11.91 12.05 6.01 5.34 
9.53 9.58 4.76 4.17 

7.27 7.27 3.60 3.09 
52.12 51.97 25.95 21.84 
35.16 34.99 17.80 14.63 
22.53 22.42 11.51 9.32 
13.86 13.68 7.09 5.66 

8.21 8.04 4.15 3.29 

4.74 4.57 2.35 1.85 

26.66 25.42 12.87 10.06 

14.92 13.96 6.84 5.34 

8.28 7.56 3.56 2.78 
4.56 4.06 1.83 1.42 

24.64 21.39 9.19 7.12 

13.16 11.16 4.54 3.48 

6.95 5.75 2.18 1.67 

36.42 29.46 10.28 7.85 

18.95 15.04 4.80 3.64 

9.78 7.62 2.19 1.67 

50.15 38.59 9.86 7.52 
25.50 19.37 4.43 3.37 
12.94 9.71 1.97 1.50 
65.41 48.57 8.70 6.72 
33.06 24.31 3.84 2.99 
16.58 12.09 1.71 1.33 

84.16 60.73 7.67 6.08 



7.34 7.43 
6.01 6.08 
4.74 4.77 

35.54 35.52 
25.33 25.16 
17.23 17.02 

11.26 11.04 
7.07 6.87 
4.29 4.11 

25.21 23.81 
14.46 13.45 
8.15 7.43 
4.53 4.02 

24.55 21.38 
13.16 11.17 

6.95 5.76 
36.45 29.54 
18.99 15.06 

9.80 7.64 
50.20 38.63 
25.50 19.34 
12.93 9.70 
65.45 48.56 
33.02 24.29 
16.62 12.08 

84.27 60.79 



3.56 3.51 

2.90 2.81 

2.26 2.14 

16.86 15.69 

12.03 10.97 

8.21 7.30 

5.35 4.65 

3.33 2.85 

2.00 1.67 

11.52 9.44 

6.42 5.15 
3.45 2.73 
1.80 1.41 

9.17 7.09 
4.53 3.47 

2.18 1.67 
10.32 7.85 

4.79 3.64 

2.20 1.67 

9.87 7.53 

4.43 3.37 
1.97 1.50 
8.70 6.70 
3.84 3.00 
1.70 1.34 
7.66 6.06 
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TABLE II. Neutrino flux (m-^sec-^sr-^GeV-^) for 0.9 > cos6l^ > 0.8 



(GeV) 



Kamioka 



Sudbury 



Gran Sasso 



Up Vp 



Norm 



4.33 4.33 
3.57 3.57 
2.84 2.84 
21.65 21.65 
15.80 15.80 
11.09 11.09 
7.53 7.53 
4.94 4.94 
3.13 3.13 
19.26 19.26 
11.52 11.52 
6.68 6.68 
3.78 3.78 
20.80 20.80 

11.19 11.19 
5.88 5.88 

30.44 30.44 
15.63 15.63 
7.94 7.94 

40.20 40.20 
20.11 20.11 
10.08 10.08 
50.50 50.50 
25.24 25.24 
12.62 12.62 
63.49 63.49 



2.06 2.11 
1.69 1.71 

1.33 1.33 
10.20 9.94 

7.47 7.17 

5.27 4.95 
3.58 3.30 

2.34 2.12 

1.48 1.31 

9.03 7.89 
5.33 4.56 

3.04 2.54 
1.68 1.38 
8.94 7.20 
4.63 3.65 
2.30 1.79 

11.14 8.58 

5.28 4.00 
2.43 1.84 

11.05 8.38 
4.98 3.77 
2.21 1.68 
9.77 7.51 
4.30 3.34 
1.91 1.49 
8.60 6.79 



12.25 12.38 
9.71 9.76 
7.37 7.38 
52.72 52.62 
35.57 35.44 
22.84 22.72 
14.09 13.95 
8.36 8.22 
4.82 4.69 
27.38 26.25 
15.29 14.43 
8.48 7.83 
4.68 4.20 
25.31 22.26 
13.49 11.60 
7.14 5.99 
37.40 30.66 
19.46 15.68 
10.03 7.93 
51.52 40.26 
26.17 20.12 
13.28 10.07 
67.24 50.40 
34.01 25.22 
17.11 12.55 
87.22 63.37 



6.17 5.50 

4.87 4.27 

3.68 3.17 

26.52 22.35 

18.22 15.02 

11.84 9.59 

7.32 5.87 

4.33 3.43 
2.47 1.94 

13.54 10.64 

7.27 5.68 

3.82 2.97 

1.97 1.54 

9.98 7.73 

4.95 3.81 
2.40 1.83 

11.36 8.65 

5.31 4.03 

2.44 1.84 

11.04 8.43 

4.96 3.76 
2.21 1.68 
9.77 7.48 
4.33 3.34 
1.90 1.50 
8.56 6.77 



7.50 7.61 
6.12 6.18 
4.81 4.83 
35.98 35.94 
25.60 25.49 
17.48 17.27 
11.41 11.22 
7.19 7.00 
4.38 4.21 
25.78 24.48 
14.82 13.87 
8.36 7.69 
4.63 4.17 
25.26 22.14 
13.48 11.60 
7.15 6.00 
37.50 30.73 
19.47 15.70 
10.03 7.97 
51.69 40.13 
26.21 20.10 
13.30 10.09 
67.29 50.42 
34.01 25.22 
17.13 12.60 
87.19 63.49 



3.65 3.59 

2.95 2.86 

2.30 2.19 

17.24 15.99 

12.32 11.24 

8.45 7.51 

5.54 4.82 

3.48 2.96 

2.10 1.75 

12.18 9.99 

6.82 5.49 

3.69 2.92 

1.94 1.52 

9.91 7.68 

4.94 3.79 

2.40 1.83 

11.38 8.67 

5.33 4.01 

2.45 1.84 

11.17 8.45 

4.97 3.76 

2.21 1.69 

9.74 7.50 

4.32 3.34 

1.91 1.49 

8.58 6.77 
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TABLE III. Neutrino flux (m-^sec-^sr-^GeV-^) for 0.8 > cos6'^ > 0.7 



(GeV) 



Kamioka 



Vp Up 



Sudbury 



Gran Sasso 



Up Up 



Norm 



4.46 4.46 2.12 2.17 

3.63 3.63 1.73 1.74 

2.88 2.88 1.36 1.35 

2.18 2.18 1.04 1.01 

15.92 15.92 7.63 7.30 

11.19 11.19 5.41 5.07 
7.60 7.60 3.68 3.40 
5.00 5.00 2.42 2.20 
3.18 3.18 1.54 1.37 

19.55 19.55 9.42 8.27 

11.76 11.76 5.62 4.81 

6.85 6.85 3.24 2.72 

3.90 3.90 1.80 1.48 

21.52 21.52 9.67 7.81 

11.62 11.62 5.03 3.98 

6.12 6.12 2.53 1.97 

31.87 31.87 12.36 9.49 

16.31 16.31 5.85 4.49 

8.30 8.30 2.73 2.08 

41.98 41.98 12.58 9.60 

21.02 21.02 5.63 4.27 

10.54 10.54 2.51 1.91 

52.84 52.84 11.19 8.51 

26.43 26.43 4.91 3.80 

13.20 13.20 2.18 1.70 
66.54 66.54 9.75 7.67 



12.68 12.82 6.38 5.70 

9.97 10.02 5.01 4.40 

7.50 7.51 3.77 3.24 

5.35 5.33 2.71 2.29 
36.02 35.91 18.63 15.42 
23.15 23.10 12.18 9.89 
14.28 14.19 7.57 6.06 

8.51 8.39 4.49 3.58 
4.92 4.82 2.58 2.04 

27.97 27.00 14.33 11.23 

15.64 14.88 7.74 6.06 

8.69 8.12 4.09 3.19 

4.81 4.36 2.13 1.66 

26.07 23.12 10.89 8.40 

13.91 12.09 5.41 4.19 

7.36 6.26 2.65 2.02 

38.67 32.12 12.65 9.61 
20.02 16.39 5.95 4.49 
10.35 8.33 2.74 2.07 

53.19 42.10 12.53 9.40 
27.00 21.02 5.65 4.26 
13.72 10.54 2.51 1.92 
69.61 52.73 11.17 8.56 

35.20 26.35 4.91 3.81 
17.76 13.17 2.18 1.69 

90.68 66.58 9.78 7.69 



7.74 7.84 
6.25 6.31 
4.88 4.90 
3.64 3.64 
25.92 25.74 
17.65 17.48 
11.59 11.41 
7.31 7.14 
4.46 4.31 
26.37 25.16 
15.18 14.31 
8.58 7.95 
4.77 4.33 
25.95 23.04 
13.92 12.11 
7.37 6.27 
38.57 32.19 
20.04 16.41 
10.36 8.31 
53.22 41.93 
27.01 21.02 
13.72 10.54 
69.59 52.87 
35.18 26.41 
17.79 13.19 
90.69 66.55 



3.75 3.70 
3.03 2.94 
2.35 2.24 

1.76 1.64 
12.63 11.50 

8.69 7.72 

5.72 4.98 

3.62 3.08 

2.20 1.84 

12.85 10.54 

7.25 5.83 

3.95 3.14 

2.10 1.64 

10.80 8.41 

5.41 4.18 

2.65 2.02 

12.67 9.62 

5.99 4.51 

2.75 2.09 

12.40 9.43 

5.65 4.27 

2.52 1.92 

11.17 8.52 

4.92 3.80 

2.18 1.70 

9.75 7.69 
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TABLE IV. Neutrino flux (m-^sec-^sr-^GeV-M for 0.7 > cos6l^ > 0.6 



(GeV) 



Kamioka 



Vp Up 



Sudbury 



Gran Sasso 



Up Up 



Norm 



4.63 4.63 2.20 2.24 

3.74 3.74 1.78 1.79 

2.94 2.94 1.40 1.38 

2.22 2.22 1.06 1.03 

16.09 16.09 7.82 7.43 

11.28 11.28 5.54 5.17 
7.69 7.69 3.79 3.49 
5.05 5.05 2.51 2.27 
3.22 3.22 1.61 1.42 

19.92 19.92 9.94 8.68 

12.01 12.01 5.97 5.08 

7.05 7.05 3.44 2.89 

4.02 4.02 1.93 1.59 

22.29 22.29 10.52 8.43 
12.07 12.07 5.52 4.37 

6.39 6.39 2.79 2.18 
3.33 3.33 1.38 1.06 

17.14 17.14 6.60 5.03 

8.72 8.72 3.11 2.33 

4.40 4.40 1.43 1.08 
22.13 22.13 6.49 4.90 
11.09 11.09 2.91 2.21 
55.59 55.59 12.94 9.83 
27.78 27.78 5.73 4.41 
13.91 13.91 2.54 1.96 
70.32 70.32 11.32 8.87 



13.31 13.43 6.67 5.99 

10.31 10.36 5.20 4.58 

7.70 7.70 3.88 3.35 

5.45 5.44 2.79 2.36 

36.61 36.53 19.16 15.84 

23.48 23.47 12.48 10.20 

14.51 14.45 7.80 6.28 

8.66 8.58 4.68 3.72 

5.03 4.93 2.70 2.13 

28.65 27.74 15.11 11.89 

16.10 15.38 8.24 6.46 

8.96 8.42 4.40 3.43 

4.95 4.55 2.31 1.79 

26.85 24.20 11.89 9.20 

14.36 12.70 6.00 4.60 

7.61 6.58 2.95 2.25 

3.98 3.38 1.42 1.08 

20.71 17.27 6.71 5.10 

10.71 8.71 3.12 2.38 

5.49 4.42 1.43 1.09 

27.97 22.12 6.49 4.91 

14.22 11.09 2.91 2.21 

72.17 55.62 12.98 9.88 

36.65 27.87 5.74 4.39 

18.50 13.92 2.53 1.97 

94.48 70.25 11.34 8.89 



8.05 8.17 

6.44 6.50 

4.98 5.02 
3.70 3.69 

26.20 26.11 
17.87 17.71 
11.73 11.55 

7.42 7.25 

4.54 4.39 

26.89 25.85 

15.57 14.77 

8.82 8.24 

4.92 4.51 

26.84 24.18 

14.34 12.66 

7.62 6.58 

3.99 3.38 
20.77 17.24 
10.69 8.73 

5.49 4.41 

27.96 22.10 

14.21 11.09 
72.11 55.59 
36.61 27.82 
18.49 13.90 
94.53 70.30 



3.90 3.85 
3.12 3.03 

2.41 2.29 
1.80 1.68 

12.95 11.77 

8.92 7.93 

5.89 5.15 

3.75 3.20 

2.29 1.92 

13.55 11.09 

7.72 6.22 

4.26 3.38 

2.27 1.78 
11.80 9.13 

5.96 4.58 

2.95 2.25 

1.42 1.08 

6.71 5.07 
3.14 2.35 

1.43 1.08 
6.48 4.90 

2.91 2.21 
12.92 9.86 

5.72 4.39 
2.52 1.95 

11.32 8.90 
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TABLE V. Neutrino flux (m-^sec-^sr-^GeV-^) for 0.6 > cos6l^ > 0.5 







Kamioka 








Sudbury 








Gran Sasso 








(GeV) 
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TABLE VI. Neutrino flux (m-^sec-^sr-^GeV-M for 0.5 > cos6l^ > 0.4 



(GeV) 



Kamioka 



Vp Up 



Sudbury 



Gran Sasso 



Up Up 



Norm 



5.28 5.28 2.49 2.52 

4.15 4.15 1.97 1.97 

3.17 3.17 1.53 1.48 

2.35 2.35 1.14 1.09 

16.79 16.79 8.31 7.81 

11.64 11.64 5.86 5.41 

7.87 7.87 4.01 3.65 

5.17 5.17 2.68 2.39 

3.31 3.31 1.73 1.53 

20.67 20.67 10.91 9.47 

12.57 12.57 6.65 5.66 

7.44 7.44 3.93 3.28 

4.31 4.31 2.25 1.86 

2.43 2.43 1.25 1.01 

13.20 13.20 6.70 5.35 

7.09 7.09 3.51 2.73 

3.74 3.74 1.76 1.36 

19.32 19.32 8.68 6.61 

9.86 9.86 4.17 3.17 

5.01 5.01 1.97 1.48 

25.18 25.18 9.07 6.82 

12.64 12.64 4.13 3.11 

6.33 6.33 1.86 1.41 

31.70 31.70 8.36 6.34 

15.87 15.87 3.69 2.84 

8.04 8.04 1.66 1.28 



15.55 15.71 7.75 6.99 

11.65 11.72 5.89 5.22 

8.45 8.46 4.30 3.73 

5.86 5.84 3.03 2.57 

38.68 38.57 20.57 17.12 

24.52 24.50 13.39 10.96 
15.06 15.04 8.36 6.75 

9.02 8.94 5.04 4.05 

5.25 5.19 2.96 2.35 

30.11 29.49 16.85 13.31 

17.01 16.54 9.36 7.35 

9.57 9.15 5.11 4.00 

5.31 5.01 2.74 2.13 

2.91 2.69 1.44 1.13 

15.58 14.20 7.45 5.76 

8.28 7.41 3.74 2.87 

4.34 3.83 1.84 1.40 

22.62 19.58 8.93 6.69 

11.65 9.97 4.21 3.17 

5.99 5.02 1.97 1.48 

30.51 25.21 9.10 6.84 

15.53 12.62 4.12 3.11 
7.90 6.32 1.86 1.41 

40.21 31.73 8.33 6.34 

20.42 15.83 3.69 2.84 

10.45 8.04 1.67 1.28 



9.26 9.39 
7.18 7.25 
5.42 5.44 
3.93 3.92 
27.36 27.26 
18.40 18.27 
12.03 11.87 
7.62 7.47 
4.68 4.55 
27.92 27.09 
16.29 15.66 

9.34 8.89 
5.24 4.92 
2.89 2.67 

15.54 14.19 
8.29 7.42 
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22.58 19.58 
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10.47 8.05 
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9.41 8.39 

6.22 5.46 
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1.43 1.11 

7.39 5.71 

3.74 2.85 

1.84 1.40 

8.90 6.73 

4.24 3.19 

1.98 1.48 

9.09 6.83 

4.13 3.12 

1.87 1.41 

8.33 6.35 

3.71 2.86 

1.66 1.29 
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TABLE VII. Neutrino flux (m-^sec-^sr-^GeV-^) for 0.4 > cos6'2 > 0.3 







Kamioka 






Sudbury 






Gran Sasso 






(GeV) 
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o.ou 


3.36 


1.78 


1 i^fi 


O.OiJ 


5.32 


3.08 


9 4fi 
z .^u 


4 79 


4.61 


2.54 


9 1 fi 

Z . lU 


1 n^ 


VJ. 1 C/'4:0 


90 SQ 


20.89 


11.26 


Q 71 


'i^ nn 

Ol.UU 


30.36 


17.76 


1 4 n4 


98 9n 


27.34 


15.38 


1 9 77 

IZ. / 1 




1 nnnn 


1 9 7^ 


12.75 


6.95 


o.oy 


1 7 ^'^ 

1 / .(JO 


17.12 


9.92 


7 8^ 


1 fi '^'^ 

lU.OO 


15.95 


9.03 


7 "^fi 
/ .ou 




1 9^RQ 


7 fin 


7.60 


4.15 


47 


Q 8S 
y .oo 


9.51 


5.51 




y .ou 


9.12 


5.15 


4 1 fi 

^. lU 


ml 




A A^ 


4.45 


2.40 


r .yo 


'i9 


5.23 


2.98 


9 "^4 
z . o^ 


41 


5.11 


2.86 


9 9Q 
z . z y 


ml 




9 ^1 


2.51 


1.36 


1 1 n 

1 . lU 


n'^ 

O.UO 


2.84 


1.59 


1 94 


9 QQ 
z.yy 


2.80 


1.56 


1 99 
1 .zz 


ml 


9 "^1 1 Q 


lO.OO 


13.88 


7.44 


u.yu 


1 fi "^fi 

lU.OU 


15.18 


8.34 


fi 4fi 
u.^u 


1 fi 9*^ 


15.00 


8.25 


fi "^Q 
u.oy 


mo 




7 ^9 


7.52 


3.93 


ns 

o.uo 


8 7fi 
o. / u 


7.97 


4.29 


97 

O.Z 1 


8 fi7 
o.u / 


7.92 


4.25 


'\ 9fi 
o.zu 


mo 


Q^l 1 


zL nn 


4.00 


2.03 




4 '\Q 


4.15 


2.14 


1 fi'^ 

l.UO 


4 1^8 


4.13 


2.13 


1 fil 


mo 


01 1 Q 


9n 


20.85 


10.09 


7 74 


9*^ Q8 


21.39 


10.48 


I .yo 


9"^ Q4 


21.27 


10.41 


7 Q4 
/ .y^ 


10-1 




1 n 7n 


10.70 


4.97 


74 




10.86 


5.08 


89 


1 9 "^"^ 


10.82 


5.04 


% 78 
o. / o 


10-1 






5.45 


2.38 


1 7Q 
X . 1 y 


fi "^Q 
u.oy 


5.53 


2.41 


1 81 
x.ox 


fi '^fi 
u.ou 


5.49 


2.39 


1 78 
X . / o 


10-1 


1 nnn 

xu.uuu 


97 ^^n 


Z ^.50 


11 1 o 

11.12 


8 "^fi 
o.ou 


■^9 A'\ 


Zi.ob 


11 IV 

ii.i 1 


8 "^7 
o.o / 


■^9 41 


2 l.bv 


11 1/1 
ii.i4 


8 "^fi 
o.ou 


1 n~2 


12.589 


13.80 


13.80 


5.13 


3.86 


16.53 


13.84 


5.13 


3.85 


16.49 


13.80 


5.12 


3.86 


10-2 


15.849 


6.91 


6.91 


2.34 


1.76 


8.43 


6.93 


2.34 


1.76 


8.41 


6.91 


2.34 


1.76 


10-2 


19.953 


34.65 


34.65 


10.60 


8.02 


42.83 


34.68 


10.55 


8.01 


42.70 


34.59 


10.49 


8.04 


10-3 


25.119 


17.35 


17.35 


4.72 


3.60 


21.72 


17.40 


4.75 


3.64 


21.73 


17.37 


4.73 


3.62 


10-3 


31.623 


8.80 


8.80 


2.14 


1.64 


11.21 


8.81 


2.14 


1.65 


11.20 


8.81 


2.14 


1.65 


10-3 
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TABLE VIII. Neutrino flux (m-^sec^^sr-^GeV^^) for 0.3 > cos6l^ > 0.2 



(GeV) 



Kamioka 



Sudbury 



Gran Sasso 



Up Vp 



Norm 



0.1000 


7.15 


7.15 


0.1259 


5.36 


5.36 


0.1585 


3.95 


3.95 


0.1995 


2.80 


2.80 


0.2512 


19.25 


19.25 


0.3162 


12.89 


12.89 


0.3981 


8.48 


8.48 


0.5012 


5.44 


5.44 


0.6310 


3.44 


3.44 


0.7943 


2.13 


2.13 


1.0000 


12.99 


12.99 


1.2589 


7.75 


7.75 


1.5849 


4.56 


4.56 


1.9953 


2.59 


2.59 


2.5119 


14.54 


14.54 


3.1623 


8.00 


8.00 


3.9811 


4.29 


4.29 


5.0119 


22.71 


22.71 


6.3096 


11.87 


11.87 


7.9433 


6.07 


6.07 


10.000 


3.07 


3.07 


12.589 


15.50 


15.50 


15.849 


7.78 


7.78 


19.953 


3.90 


3.90 


25.119 


19.54 


19.54 


31.623 


9.91 


9.91 



3.32 3.36 
2.55 2.52 
1.90 1.83 
1.38 1.30 
9.74 8.95 
6.68 6.04 
4.47 3.93 
2.93 2.54 
1.89 1.61 
1.18 1.00 
7.25 6.09 
4.37 3.62 
2.57 2.10 
1.47 1.18 
8.20 6.54 
4.41 3.50 
2.34 1.82 
11.98 9.23 
6.03 4.57 
2.95 2.23 
1.41 1.06 
6.63 4.98 
3.08 2.32 
1.41 1.07 
6.41 4.88 
2.93 2.25 



21.70 21.92 
15.49 15.62 
10.79 10.79 
7.16 7.13 
45.31 45.26 
27.75 27.67 
16.56 16.50 
9.69 9.64 
5.60 5.51 

3.20 3.13 
18.10 17.74 
10.25 9.90 

5.73 5.48 
3.16 3.01 
17.12 16.13 

9.21 8.61 
4.88 4.51 

25.63 23.36 

13.33 11.95 
6.84 6.08 
3.50 3.08 

17.86 15.49 
9.12 7.78 
4.64 3.90 

23.66 19.52 



10.74 9.74 

7.86 7.01 

5.51 4.81 

3.75 3.21 

24.59 20.55 

15.46 12.78 
9.42 7.69 
5.59 4.51 
3.25 2.62 
1.86 1.49 

10.47 8.34 
5.85 4.63 
3.23 2.53 
1.75 1.36 
9.34 7.23 
4.89 3.74 
2.50 1.90 

12.51 9.45 

6.15 4.63 

2.98 2.23 

1.41 1.06 

6.63 4.97 

3.07 2.30 

1.41 1.07 

6.39 4.90 



12.20 9.91 2.93 2.25 



12.65 12.85 
9.37 9.46 
6.79 6.81 
4.72 4.70 

31.50 31.30 

20.40 20.22 

12.91 12.70 
7.96 7.81 
4.81 4.70 
2.86 2.78 

16.77 16.17 
9.66 9.29 
5.51 5.25 
3.08 2.91 

16.89 15.84 
9.12 8.47 
4.84 4.46 

25.44 23.26 

13.17 11.94 
6.88 6.11 
3.50 3.08 

17.88 15.49 
9.12 7.78 
4.64 3.90 

23.61 19.54 

12.18 9.91 



6.00 5.98 

4.52 4.44 

3.30 3.15 

2.34 2.18 

16.08 14.67 

10.68 9.50 

6.87 6.02 

4.32 3.71 
2.66 2.24 
1.59 1.33 

9.33 7.68 

5.39 4.39 
3.03 2.44 
1.68 1.33 
9.06 7.14 
4.79 3.72 
2.47 1.88 

12.42 9.36 

6.14 4.62 

2.99 2.23 

1.41 1.06 

6.63 4.97 

3.08 2.32 

1.41 1.06 

6.40 4.88 
2.92 2.25 
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TABLE IX. Neutrino flux (m-^sec-^sr-^GeV-M for 0.2 > cos6l^ > 0.1 



(GeV) 



Kamioka 



Sudbury 



Gran Sasso 



Norm 



9.80 
7.14 
5.08 
3.51 
2.34 

15.13 
9.58 
6.03 
3.71 
2.25 

13.61 
8.08 
4.70 
2.70 

15.27 
8.47 
4.60 
2.49 

13.08 
6.88 
3.49 

17.81 
9.03 
4.57 

22.88 

11.61 



9.80 4.53 
7.14 3.39 
5.08 2.45 
3.51 1.73 
2.34 1.19 

15.13 7.94 
9.58 5.20 
6.03 3.31 
3.71 2.09 
2.25 1.28 

13.61 7.86 
8.08 4.67 
4.70 2.76 
2.70 1.58 

15.27 8.92 
8.47 4.94 
4.60 2.67 
2.49 1.42 

13.08 7.26 
6.88 3.70 
3.49 1.82 

17.81 8.85 
9.03 4.22 
4.57 1.99 

22.88 9.32 

11.61 4.35 



4.53 
3.32 
2.35 
1.62 
1.09 
7.04 
4.49 
2.82 
1.75 
1.06 
6.40 
3.77 
2.20 
1.27 
7.02 
3.88 
2.07 
1.08 
5.56 
2.80 
1.37 
6.64 
3.18 
1.51 
7.08 
3.31 



29.97 
20.86 
14.07 
9.06 
5.57 
32.95 
19.04 
10.85 
6.09 
3.40 
19.05 
10.68 
5.98 
3.31 
18.01 
9.78 
5.22 
2.74 
14.46 
7.49 
3.85 
19.76 
10.12 
5.16 
26.26 
13.63 



30 


35 


14.81 


20 


98 


10.57 


14 


11 


7.23 


9 


04 


4.78 


5 


56 


3.05 


32 


93 


18.66 


19 


02 


11.07 


10 


77 


6.35 


6 


01 


3.62 


3 


34 


2.02 


18 


61 


11.28 


10 


32 


6.28 


5 


72 


3.46 


3 


15 


1.92 


17 


12 


10.40 


9 


24 


5.55 


4 


91 


2.91 


2 


59 


1.49 



13.46 7.61 

6.95 3.76 

3.52 1.84 

17.85 8.88 

9.03 4.24 

4.54 1.99 

22.87 9.24 

11.67 4.35 



13.52 
9.47 
6.34 
4.09 
2.57 

15.46 
8.99 
5.15 
2.88 
1.62 
8.98 
4.95 
2.75 
1.50 
8.04 
4.25 
2.21 
1.13 
5.67 
2.85 
1.38 
6.66 
3.19 
1.51 
7.05 
3.32 



17.36 17.68 
12.52 12.61 
8.78 8.82 
5.92 5.92 

3.83 3.79 
23.90 23.76 
14.61 14.46 

8.75 8.64 
5.14 5.04 
3.00 2.93 

17.28 16.79 
9.96 9.55 
5.65 5.41 
3.17 3.00 

17.46 16.55 
9.53 8.99 
5.13 4.81 
2.73 2.55 

14.23 13.32 
7.46 6.87 

3.84 3.51 
19.76 17.87 
10.12 9.04 

5.17 4.55 
26.41 22.86 
13.64 11.64 



8.17 8.18 

6.02 5.92 

4.27 4.09 

2.95 2.75 

1.97 1.79 

12.75 11.32 

7.99 6.91 

4.86 4.14 

2.92 2.45 

1.72 1.42 

9.97 8.15 

5.68 4.61 

3.21 2.59 

1.79 1.44 

9.79 7.85 

5.32 4.17 

2.82 2.18 

1.46 1.12 

7.43 5.59 

3.72 2.84 

1.84 1.38 

8.89 6.66 

4.23 3.17 

1.99 1.50 

9.30 7.04 

4.34 3.33 
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TABLE X. Neutrino flux (m-^sec-^sr-^GeV-^) for 0.1 > cos6l^ > 0.0 







Kamioka 






Sudbury 








Gran Sasso 






(GeV) 


























Norm 


n 1 nnn 


1 Q4 


15.94 


7.20 


7 '^9 
/ .o^ 


48 (T^ 
40. uo 


48.70 


23.42 


91 7n 

^ 1 . 1 u 


28 


03 


28.57 


13.04 13.14 


1 

lU 


n 1 9^Q 


1 1 91 


11.21 


5.25 


^ 9n 


■^9 "^Q 
oz.oy 


32.64 


16.36 


1 4 78 
14. / 


19 


52 


19.71 


9.35 


9.21 


1 n3 

lU 


1 


7 7n 

/ . / u 


7.70 


3.68 


o.oo 


91 1 Q 
zi . ly 


21.31 


10.82 


Q fi4 


13 


23 


13.31 


6.45 


6.20 


1 n3 

lU 




tj. lO 


5.13 


2.51 




1 ^ 1 Q 
10. ly 


13.19 


6.97 


09 


8 


61 


8.56 


4.30 


4.01 


1 n3 

lU 






3.27 


1.67 


l.ol 


7 7Q 
1 . 1 y 


7.81 


4.31 


o.uo 


5 


34 


5.31 


2.78 


2.51 


1 n3 

lU 




90 1 1 

. 1 1 


20.11 


10.78 


Q 48 


44 4fi 

4:4:. 4:U 


44.36 


25.43 


91 9n 


31 


96 


31.86 


17.31 


15.26 


102 

lU 




1 9 "^n 


12.30 


6.84 


^ 74 

O. 1 4 


9zL fifi 


24.57 


14.53 


1 1 8^ 
1 1.00 


18 


69 


18.51 


10.44 


8.96 


1 n^ 

lU 


n "^01 9 


7 4zL 


7.44 


4.17 


0.4U 


10. 4:1 


13.43 


8.11 


4Q 


10 


78 


10.61 


6.08 


5.19 


1 n^ 

lU 


W.UOIU 




4.40 


2.53 


9 ns 


7 9fi 


7.15 


4.38 


o.oo 


6 


08 


5.92 


3.57 


2.93 


1 n^ 

lU 




9 ^8 


2.58 


1.52 


1 91 
i.zi 


o.yo 


3.87 


2.41 


1 on 

1 . yu 


3 


40 


3.32 


2.03 


1.65 


1 n^ 

lU 


1 nnnn 


1 1 7 
itj. 1 / 


15.17 


9.00 


7 1 n 

/ . lU 


91 9"^ 


20.81 


12.85 


1 9*^ 


19 


06 


18.64 


11.33 


9.18 


101 

lU 


1 9^RQ 


o.oo 


8.83 


5.25 


A 1 1 
4.11 


1 1 (^7 


11.27 


6.97 


0.00 


10 


71 


10.38 


6.31 


5.09 


101 

lU 




O. lU 


5.10 


3.08 


9 '^7 

^.O 1 


41 

U .4:X 


6.17 


3.77 


9 QQ 
^.yy 


5 


98 


5.76 


3.52 


2.78 


101 

lU 




9 87 


2.87 


1.75 


1 .00 


^ ^9 


3.38 


2.08 


1 .uo 


3 


32 


3.19 


1.95 


1.54 


101 

lU 


9 "^1 1 Q 


IRIS 

lU. lO 


16.18 


9.73 


7 zLQ 


1 Q n'^ 

ly .uo 


18.31 


11.29 


8 8*^ 
0.00 


18 


47 


17.45 


10.66 


8.43 


100 

lU 




O.c/O 


8.93 


5.45 


4 1 Q 

4. 1 c? 




9.78 


6.04 


4 fi8 

4:.UO 


9 


94 


9.47 


5.84 


4.54 


10° 

lU 


Q^l 1 


zL Q1 


4.91 


2.96 


9 97 


0.01 


5.22 


3.22 


9 47 


5 


41 


5.13 


3.07 


2.43 


10° 

lU 


01 1 Q 


9 RA 


2.64 


1.59 


1 9n 


9 Q9 


2.79 


1.69 


1 9Q 
i.zy 


2 


88 


2.75 


1.66 


1.28 


10° 

lU 




1 zL 9fi 


14.26 


8.31 


u.oo 


10.00 


14.48 


8.83 


fi9 


15 


38 


14.55 


8.73 


6.69 


10-1 

lU 




7 47 


7.47 


4.34 


o.ou 


o.uy 


7.65 


4.48 


4"^ 
0.40 


8 


(TV 

07 


7.64 


4.56 


3.47 


10-1 


1 nnn 

xu.uuu 


O.c/O 


O HQ 

o.yo 


z.zo 


1 f!Q 


A 91 


4.00 


z.oO 


1 79 
1 . / z 


4 


20 


o.9d 


2.28 


1.72 


10-1 

lU 


12.589 


20.27 


20.27 


11.34 


8.50 


21.81 


20.55 


11.52 


8.66 


21 


82 


20.49 


11.50 


8.67 


10-2 


15.849 


10.48 


10.48 


5.66 


4.25 


11.31 


10.57 


5.69 


4.27 


11 


33 


10.58 


5.71 


4.26 


10-2 


19.953 


5.39 


5.39 


2.78 


2.08 


5.83 


5.36 


2.78 


2.11 


5 


83 


5.40 


2.79 


2.11 


10-2 


25.119 


2.73 


2.73 


1.35 


1.02 


3.01 


2.71 


1.36 


1.03 


2 


99 


2.73 


1.36 


1.02 


10-2 


31.623 


14.02 


14.02 


6.52 


4.98 


15.65 


14.08 


6.54 


5.01 


15 


62 


14.06 


6.57 


5.02 


10-3 
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TABLE XI. Neutrino flux (m-^sec-^sr-^GeV-^) for 0.0 > cos6l^ > -0.1 







Kamioka 






Sudbury 








Gran Sasso 






(GeV) 










H' 
















Norm 


n 1 nnn 


17 11 


17.11 


7.76 


7 8fi 


4fi 90 


46.96 


22 AA 


9n 


28 


33 


28.90 


13.36 13.18 


1 n^ 


n 1 9^Q 


1 1 Q7 


11.97 


5.60 


O.OiJ 


■^1 94 


31.62 


15.76 


1 4 "^8 


19 


72 


19.87 


9.46 


9.26 


1 n3 


1 


S 1 8 


8.18 


3.92 


78 


9n fin 


20.68 


10.48 


Q 4n 

y .^u 


13 


29 


13.37 


6.53 


6.18 


1 n3 




o.ou 


5.36 


2.65 


9 4Q 


1 9 Q9 


12.89 


6.80 


^ Q9 


8 


58 


8.55 


4.33 


3.98 


1 n3 




o.oo 


3.38 


1.73 




7 fi8 


7.67 


4.21 


o.oy 


5 


30 


5.28 


2.78 


2.50 


1 n3 




^u.uo 


20.65 


11.17 


Q 8^ 


44 Dfi 


44.13 


25.15 


91 n'^ 


31 


55 


31.53 


17.27 15.00 


in2 

-LU 




1 9 


12.58 


6.95 


^ Q1 


94 fifi 


24.54 


14.53 


1 1 81 


18 


31 


18.35 


10.35 


8.81 


1 n^ 


n "^01 9 


7 49 


7.42 


4.22 


'i 1^9 


1 ^1 


13.46 


8.12 


fi 1^8 


10 


49 


10.43 


6.03 


5.09 


1 n^ 


W.UOIU 




4.40 


2.54 


9 nfi 


7 "^8 


7.28 


4.44 


^7 

O.O 1 


5 


97 


5.85 


3.49 


2.87 


1 n^ 






2.56 


1.51 


1 9*^ 


4 nn 


3.93 


2.43 


1 Q1 
i.yr 


3 


36 


3.29 


2.01 


1.62 


1 n^ 


1 nnnn 


1 4 Qfi 


14.96 


8.92 


7 


91 ^n 


21.15 


13.02 


1 n "^Q 
lu.oy 


18 


83 


18.40 


11.17 


9.00 




1 9^RQ 


S fi7 


8.67 


5.23 


4 (T^ 


1 1 78 


11.34 


7.12 


o.uo 


10 


62 


10.33 


6.28 


5.07 






^ 04 


5.04 


3.02 


9 '^8 


fi 'in 


6.22 


3.83 


'\ n4 


5 


90 


5.76 


3.51 


2.77 


ml 

-LU 




9 S^^ 


2.85 


1.72 




1^4 


3.40 


2.08 


1 fii^ 


3 


32 


3.19 


1.96 


1.54 


ml 

lU 


9 "^1 1 Q 


1 nfi 

lU.UU 


16.06 


9.56 


7 44 


10 91^ 
ly.Ztj 


18.41 


11.19 


8 Q1 


18 


41 


17.46 


10.75 


8.42 


mo 

lU 




S Q9 


8.92 


5.41 


4 14 


1 n '^fi 


9.83 


6.11 


4 79 


10 


02 


9.49 


5.87 


4.53 


mo 

lU 


Q^l 1 


4 87 


4.87 


2.91 


9 97 


^ ^1 


5.27 


3.23 


9 47 


5 


43 


5.15 


3.10 


2.41 


mo 

lU 


01 1 Q 




2.63 


1.58 


1 91 


9 Q9 


2.81 


1.71 


1 9Q 


2 


90 


2.73 


1.64 


1.28 


mo 

lU 




1 4 nn 


14.00 


8.33 


u.ou 


1 4"^ 


14.55 


8.88 


fi fi4 


15 


31 


14.45 


8.69 


6.58 


in-i 




7 4Q 


7.49 


4.35 


41 


8 19 


7.63 


4.52 


% 44 


8 


12 


7.66 


4.45 


3.41 


10-1 


1 nnn 

xu.uuu 


O.iJO 


0.9o 


z.z4 


1 f!Q 


4 99 


0.99 


1.61 


1 79 


4 


20 


o.9o 


2.28 


1.73 


10-1 


12.589 


20.29 


20.29 


11.34 


8.48 


21.88 


20.55 


11.57 


8.62 


21 


84 


20.51 


11.49 


8.67 


10-2 


15.849 


10.48 


10.48 


5.67 


4.26 


11.30 


10.57 


5.71 


4.27 


11 


35 


10.59 


5.72 


4.25 


10-2 


19.953 


5.40 


5.40 


2.78 


2.08 


5.82 


5.38 


2.78 


2.11 


5 


83 


5.39 


2.80 


2.12 


10-2 


25.119 


2.73 


2.73 


1.36 


1.02 


3.01 


2.72 


1.35 


1.03 


2 


99 


2.73 


1.37 


1.02 


10-2 


31.623 


14.03 


14.03 


6.50 


4.97 


15.66 


14.05 


6.55 


5.00 


15 


60 


14.09 


6.58 


5.01 


10-3 
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TABLE XII. Neutrino flux (m-^sec-^sr-iGeV"^) for -0.1 > cosO^ > -0.2 



(GeV) 



Kamioka 



Sudbury 



Gran Sasso 



Norm 



12.56 12.56 
8.94 8.94 
6.24 6.24 
4.18 4.18 
2.69 2.69 

16.89 16.89 

10.41 10.41 
6.32 6.32 
3.81 3.81 
2.27 2.27 

13.47 13.47 
7.93 7.93 
4.61 4.61 
2.66 2.66 

15.07 15.07 
8.36 8.36 
4.56 4.56 
2.45 2.45 

12.98 12.98 
6.77 6.77 
3.49 3.49 

17.79 17.79 
9.01 9.01 
4.54 4.54 

22.74 22.74 



5.81 5.81 
4.27 4.18 

3.04 2.89 
2.09 1.94 
1.40 1.27 

9.05 8.03 
5.71 4.95 

3.53 3.02 
2.17 1.81 
1.31 1.09 

7.80 6.37 

4.62 3.75 
2.69 2.17 

1.54 1.25 
8.71 6.98 
4.86 3.83 

2.63 2.05 
1.39 1.08 
7.27 5.50 
3.68 2.80 

1.81 1.36 
8.85 6.63 
4.23 3.18 
2.00 1.50 



9.23 7.02 
11.60 11.60 4.35 3.30 



25.56 25.94 
17.98 18.12 
12.30 12.38 
8.07 8.06 
5.05 5.05 

30.64 30.52 
18.09 18.01 
10.47 10.39 

6.00 5.93 
3.41 3.35 

19.13 18.73 

10.76 10.47 
6.03 5.81 
3.34 3.19 

18.18 17.30 
9.79 9.28 
5.20 4.94 
2.75 2.57 

14.34 13.44 
7.49 6.95 
3.84 3.52 

19.74 17.84 

10.12 9.01 
5.17 4.54 

26.41 22.79 

13.65 11.65 



12.51 11.63 
9.03 8.25 
6.24 5.59 
4.20 3.67 
2.72 2.33 

17.09 14.31 

10.35 8.52 

6.08 4.98 

3.53 2.84 

2.01 1.61 

11.35 9.04 

6.34 5.03 
3.52 2.77 
1.93 1.50 

10.52 8.05 
5.55 4.26 

2.90 2.22 
1.49 1.13 
7.55 5.71 
3.74 2.81 
1.84 1.37 

8.91 6.64 

4.23 3.19 
2.00 1.51 

9.24 7.02 

4.35 3.31 



17.71 18.04 8.58 8.18 

12.60 12.72 6.26 5.84 

8.69 8.78 4.36 3.99 

5.77 5.80 2.97 2.65 

3.67 3.68 1.95 1.71 

22.59 22.71 12.46 10.64 

13.69 13.71 7.67 6.45 

8.16 8.13 4.61 3.86 
4.81 4.77 2.75 2.27 
2.81 2.76 1.63 1.32 

16.33 16.06 9.45 7.67 

9.47 9.19 5.44 4.40 

5.45 5.25 3.11 2.49 

3.08 2.95 1.75 1.40 

17.20 16.37 9.73 7.71 

9.45 8.93 5.28 4.11 

5.12 4.82 2.81 2.16 

2.73 2.56 1.46 1.12 

14.32 13.29 7.43 5.70 

7.49 6.83 3.74 2.81 

3.84 3.52 1.84 1.37 

19.76 17.86 8.89 6.65 

10.11 9.03 4.23 3.18 

5.17 4.56 2.00 1.51 
26.50 22.84 9.36 7.06 
13.64 11.62 4.34 3.31 
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TABLE XIII. Neutrino flux (m-^sec-isr-iGeV-i) for -0.2 > cos6l^ > -0.3 







Kamioka 






Sudbury 






Gran Sasso 






(GeV) 




























Norm 


1 nnn 




10.41 


4.95 


4 7^ 


lO.OO 


15.62 


7.58 


t .uo 


1 9 4^ 
1 ^ .^o 


12.68 


6 


14 


^ 71 

O. 1 1 


1 

lU 




/ .(JO 


7.55 


3.68 


'\ 4S 


1 1 1 Q 

1 1 . Ic/ 


11.33 


5.63 


5.15 


Q 07 
y.u / 


9.19 


4 


55 


4 1 fi 

^. lU 


1 n3 

lU 


U. l(JOO 




5.34 


2.64 


9 44 


7 0*^ 


8.00 


4.01 


3.60 


R AA 


6.49 


3 


26 




1 n3 

lU 


1 QQ"^ 

U. LlJuO 


o.uo 


3.65 


1.85 


1 fi7 


1^ 4n 


5.42 


2.79 


2.45 


4 '^Q 


4.41 


2 


26 


9 nn 


1 n3 

lU 






2.39 


1.25 


1 1 1 

1.11 


o.oo 


3.54 


1.88 


1.62 


9 88 
z.oo 


2.89 


1 


53 


1 .oo 


1 n3 

lU 




1 ^ 97 


15.27 


8.20 


7 1 ^ 
I . i(j 


99 '^8 


22.41 


12.23 10.35 


1 8 '^1 


18.37 


9 


97 


o.ou 


102 




Q f!9 


9.62 


5.24 


4 ^1 


1 S4 


13.86 


7.69 


6.45 


1 1 41 
11.^1 


11.40 


6 


32 


^ 97 
o.z / 


1 n^ 


n "^01 9 


^ Q1 


5.91 


3.29 


9 77 


8 41 


8.37 


4.75 


3.93 


fi Q7 
u.y / 


6.95 


3 


94 


9^ 


1 n^ 


U.UOIU 


O.Ul 


3.61 


2.02 


1 71 
1 . / 1 


o.uo 


4.97 


2.86 


2.33 


4 99 


4.19 


2 


38 


1 QR 
i.yu 


1 n^ 


7Q4^ 


9 17 


2.17 


1.23 


1 (T^ 


9 

z.yo 


2.91 


1.69 


1.37 


9 ^1 

Z.Ol 


2.49 


1 


43 


1 1 fi 

1 . lU 


1 n^ 


1 nnnn 


1 9 Q7 


12.97 


7.36 


U.UU 


17 18 

1 1 . lO 


16.80 


9.83 


7.90 


1 4 88 


14.61 


8 


39 


8"^ 
u.oo 


loi 




7 fifi 


7.66 


4.32 


O.fJU 


y .oo 


9.62 


5.63 


4.47 


8 74 


8.48 


4 


89 


Q4 


loi 




4 44 


4.44 


2.52 




^ fil 

O.Ul 


5.40 


3.15 


2.49 


o.uo 


4.86 


2 


79 


9 9"^ 
z . zo 


loi 


L.Vc/OO 




2.55 


1.43 


1 1 

1 . lU 


1 4 

0. 1^ 


2.98 


1.74 


1.35 


9 87 


2.73 


1 


56 


1 9^ 

l.ZO 


loi 


9 "^1 1 Q 


1 zL 9A 


14.24 


8.02 


fi 41 


1 7 07 
1 / .u / 


16.10 


9.31 


7.22 


1 ^ Q1 
lo.y 1 


15.06 


8 


62 


79 


10° 


1 R9^ 




7.84 


4.33 




Q 1 Q 
y . ly 


8.56 


4.87 


3.74 


8 78 


8.21 


4 


61 


o.ou 


10° 


QXl 1 


A 9*^ 


4.23 


2.29 


1 7S 
1 . / o 


4 88 


4.51 


2.49 


1.90 


4 79 


4.37 


2 


39 


1 8*^ 
1 .oo 


10° 


01 1 Q 


99 zL7 


22.47 


11.85 


Q 07 




23.28 


12.54 


9.47 


91^ 07 


23.04 


12 


25 


Q 9fi 
y .zu 


10-1 


u.ouyu 


1 1 fiQ 


11.69 


5.95 




1 9^ 


12.01 


6.17 


4.63 


1 1 8 

lO. lO 


11.89 


6 


10 


4 


10-1 




(\ n^ 


6.05 


2.91 


9 99 


fi 8Q 
u.oy 


6.12 


2.97 


2.25 


fi 7Q 
u. ( y 


6.06 


2 


97 


9 91 


10-1 


1 n nnn 


nfi 


o.Ub 


1 /in 
1.40 


l.UO 


o.ou 


Q no 
O.Oo 


1.42 


1.06 


'\ 4Q 


3.0/ 


1 


A 1 

4i 


1 .uu 


10-1 


12.589 


15.44 


15.44 


6.61 


4.96 


17.88 


15.50 


6.63 


4.96 


17.86 


15.44 


6 


64 


4.97 


10-2 


15.849 


7.75 


7.75 


3.07 


2.31 


9.11 


7.77 


3.07 


2.31 


9.12 


7.77 


3 


08 


2.31 


10-2 


19.953 


3.89 


3.89 


1.40 


1.06 


4.65 


3.90 


1.41 


1.07 


4.65 


3.89 


1 


41 


1.06 


10-2 


25.119 


19.46 


19.46 


6.38 


4.86 


23.63 


19.52 


6.43 


4.87 


23.65 


19.54 


6 


39 


4.87 


10-3 


31.623 


9.90 


9.90 


2.93 


2.23 


12.17 


9.90 


2.92 


2.24 


12.16 


9.88 


2 


92 


2.23 


10-3 
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TABLE XIV. Neutrino flux (m-^sec-^sr-^GeV-^) for -0.3 > cos6l^ > -0.4 







Kamioka 






Sudbury 






Gran Sasso 






(GeV) 


















H' 








Norm 


n 1 nnn 


o.ou 


8.80 


4.26 


4 ni 


1 n 71 


10.89 


5.29 


4 Q1 


Q 8Q 
y .oy 


10.05 


4.89 


4 ^4 


1 




U.fJO 


6.53 


3.21 


9 QQ 


7 QQ 


8.08 


4.00 


% fi7 


7 "^7 
1 .0 / 


7.47 


3.70 


0.00 


1 n3 




zL 79 


4.72 


2.35 


9 1^ 


'\ S9 


5.87 


2.93 


9 


(J.O / 


5.41 


2.71 


9 ZL9 


1 n3 




9Q 


3.29 


1.67 


1 


4 ns 


4.09 


2.09 


1 84 


74 


3.76 


1.92 


1 fiQ 


1 n3 




9 91 


2.21 


1.15 


1 ni 


9 74 


2.74 


1.44 


1 91^ 


9 ^1 


2.52 


1.32 


1 1 d 


1 n3 


'^1 fi9 




14.38 


7.64 


u.uu 


1 7 81 


17.86 


9.55 


8 1 8 


1 fi 98 


16.34 


8.74 


7 47 

/ .'4: 1 


102 

lU 




Q 9n 


9.20 


4.96 


4 9fi 


1 1 "if^ 


11.34 


6.18 


^ 91 


1 n "^^ 


10.38 


5.63 


4 7f^ 


1 n^ 

lU 


n "^01 9 


^ 71 


5.71 


3.12 


9 f!7 




7.01 


3.88 


% 9fi 


4^ 


6.42 


3.54 


9 Q"^ 


1 n^ 

lU 


W.UOIU 


^1 


3.51 


1.94 


1 f\A 


4 "^1 


4.26 


2.38 


1 Q7 
1 .y 1 


Q4 


3.89 


2.17 


1 7Q 


1 n^ 

lU 


VJ. 1 C/'4:0 


91 9^ 


21.25 


11.80 


Q Q1 


9^ 88 


25.41 


14.29 


1 1 7Q 
1 1 . 1 y 


9'^ 7*^ 


23.29 


13.07 


1 n 7f^ 

lU. / u 


lU 


1 nnnn 


1 9 fifi 


12.66 


7.02 


o.ou 


1 "^9 


14.90 


8.40 


fi 8^ 


1 4 HQ 


13.75 


7.69 


(\ *^n 


lU 




7 4^ 


7.45 


4.10 


41 


8 Q9 


8.62 


4.84 


QD 
0. yu 


8 94 


7.96 


4.45 


fil 

O.Ul 


ml 

lU 






4.30 


2.35 


1 Q9 


(J. LO 


4.88 


2.71 


9 1 7 


4 7?^ 


4.55 


2.50 


9 n9 


ml 

lU 




9 44 


2.44 


1.32 


1 07 


9 88 


2.71 


1.49 


1 1 7 
J — Li 


9 71 


2.55 


1.39 


1 1 1 

1.11 


ml 

lU 


9 "^1 1 Q 


1 1^9 


13.52 


7.25 


^ 7fi 


1 ^ 89 


14.68 


7.99 


fi 91 


1^0^ 
10. uo 


13.97 


7.52 


^ Q4 

tJ.y'4: 


mo 

lU 




7 "^i^ 


7.35 


3.84 


'i ni 


0.00 


7.80 


4.16 


90 


R 9*^ 


7.54 


3.96 


no 

o.uy 


mo 

lU 


Q^l 1 


O.c/O 


3.93 


1.98 


1 1^4 


4 1^4 


4.09 


2.10 


1 fin 




4.00 


2.02 


1 .ou 


10° 

lU 


01 1 Q 


9n 


20.58 


9.94 


7 fi7 


9*^ 7f\ 


21.14 


10.41 


7 8Q 
1 .oy 


9^^ A7 


20.89 


10.16 


7 74 


10-1 

lU 




1 n fi^ 


10.65 


4.91 


7^ 


1 9 "^1 


10.78 


5.07 


78 
0. i 


1 9 99 


10.73 


4.98 


*^ 7^^ 
0. / 


10-1 

lU 




^ 44 


5.44 


2.34 


1 7fi 


fi "^7 


5.46 


2.40 


1 7Q 
± . 1 y 


U.O'4: 


5.47 


2.38 


1 ?^n 

1 .ou 


10-1 

lU 


1 nnn 

xu.uuu 


97 A1 




11.06 


8 "^1 
0.01 


■^9 49 




ii.lD 


8 'W 
0.00 


^^9 '^'^ 


2 1 Ai 


11 11 
ii.ii 


R '^d 


1 n~2 

lU 


12.589 


13.79 


13.79 


5.12 


3.85 


16.52 


13.78 


5.13 


3.86 


16.50 


13.80 


5.14 


3.85 


10-2 


15.849 


6.91 


6.91 


2.33 


1.76 


8.41 


6.93 


2.34 


1.76 


8.41 


6.93 


2.34 


1.77 


10-2 


19.953 


34.66 


34.66 


10.55 


8.03 


42.75 


34.56 


10.53 


7.98 


42.78 


34.69 


10.55 


8.01 


10-^ 


25.119 


17.34 


17.34 


4.74 


3.62 


21.71 


17.35 


4.72 


3.64 


21.75 


17.33 


4.73 


3.62 


10-^ 


31.623 


8.78 


8.78 


2.13 


1.64 


11.18 


8.79 


2.14 


1.65 


11.18 


8.80 


2.14 


1.65 


10-^ 
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TABLE XV. Neutrino flux (m-^sec-^sr-^GeV-^) for -0.4 > cos 6'^ > -0.5 



(GeV) 



Kamioka 



Sudbury 



Gran Sasso 



Up Vp 



Norm 



7.89 7.89 
6.00 6.00 

4.43 4.43 
3.14 3.14 

21.51 21.51 
14.25 14.25 
9.19 9.19 
5.79 5.79 
3.56 3.56 
21.45 21.45 
12.75 12.75 

7.44 7.44 
4.25 4.25 

23.82 23.82 
13.03 13.03 
6.99 6.99 
3.68 3.68 
19.20 19.20 
9.82 9.82 

5.03 5.03 
25.16 25.16 
12.65 12.65 

6.33 6.33 
31.71 31.71 
15.85 15.85 

8.04 8.04 



3.83 3.61 

2.95 2.75 
2.20 2.01 
1.59 1.43 

11.07 9.81 
7.46 6.55 
4.89 4.23 
3.11 2.67 
1.92 1.64 

11.62 9.82 
6.87 5.75 
3.97 3.29 
2.24 1.83 

12.40 9.95 
6.59 5.29 
3.43 2.69 
1.74 1.35 
8.58 6.58 
4.13 3.15 

1.96 1.47 
9.06 6.81 
4.13 3.11 
1.86 1.41 
8.32 6.33 
3.68 2.85 
1.66 1.28 



8.34 8.46 

6.38 6.45 

4.74 4.77 

3.39 3.39 
23.24 23.29 
15.41 15.38 

9.97 9.92 
6.30 6.25 
3.88 3.83 
23.59 22.98 
14.05 13.57 
8.24 7.87 

4.75 4.47 
26.75 24.83 
14.75 13.48 

8.01 7.17 

4.26 3.76 
22.34 19.42 
11.61 9.91 

5.97 5.03 
30.53 25.21 
15.53 12.64 

7.92 6.34 
40.20 31.70 
20.45 15.90 
10.47 8.05 



4.11 3.86 

3.18 2.94 
2.37 2.16 
1.71 1.53 

12.01 10.54 

8.10 7.02 

5.30 4.54 

3.36 2.86 
2.08 1.74 

12.53 10.42 

7.37 6.07 
4.22 3.46 
2.37 1.91 

12.98 10.40 

6.92 5.44 

3.56 2.76 

1.79 1.37 

8.74 6.64 

4.19 3.15 
1.99 1.47 
9.08 6.82 
4.14 3.11 
1.86 1.41 
8.33 6.32 
3.70 2.86 
1.66 1.29 



8.52 8.64 
6.50 6.57 
4.82 4.86 
3.44 3.45 
23.42 23.48 
15.45 15.45 
9.96 9.91 
6.25 6.21 
3.85 3.78 
23.19 22.62 
13.75 13.26 
8.04 7.68 
4.64 4.36 
26.11 24.26 
14.40 13.17 
7.84 7.05 
4.20 3.70 
22.13 19.24 
11.56 9.87 
5.95 5.00 
30.47 25.16 
15.53 12.64 
7.90 6.33 
40.29 31.71 
20.44 15.91 
10.47 8.05 



4.22 3.92 

3.25 2.98 

2.42 2.18 

1.75 1.55 
12.15 10.60 

8.17 7.01 

5.32 4.51 

3.35 2.83 

2.06 1.72 

12.34 10.22 

7.22 5.95 

4.14 3.38 
2.32 1.87 

12.58 10.08 

6.74 5.30 

3.47 2.71 

1.76 1.35 
8.62 6.59 

4.15 3.15 
1.96 1.48 
9.06 6.82 
4.13 3.11 
1.86 1.41 
8.34 6.34 
3.70 2.83 
1.66 1.29 
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TABLE XVI. Neutrino flux (m-^sec-^sr-^GeV-M for -0.5 > cos6l^ > -0.6 



(GeV) 



Kamioka 



Sudbury 



Gran Sasso 



Norm 



7.89 7.89 
6.10 6.10 
4.59 4.59 
3.31 3.31 
2.28 2.28 
15.22 15.22 
9.84 9.84 
6.18 6.18 
3.77 3.77 

2.25 2.25 
13.18 13.18 

7.56 7.56 

4.26 4.26 
23.41 23.41 
12.69 12.69 

6.70 6.70 
3.50 3.50 
18.08 18.08 

9.27 9.27 
4.70 4.70 

23.43 23.43 
11.77 11.77 

5.89 5.89 
29.49 29.49 
14.75 14.75 

7.48 7.48 



3.83 3.62 
2.99 2.80 
2.26 2.08 

1.65 1.49 
1.16 1.03 
7.91 6.90 

5.16 4.46 
3.26 2.80 
2.00 1.69 
1.20 1.00 
6.95 5.75 
3.94 3.21 

2.17 1.76 
11.66 9.35 

6.13 4.81 
3.13 2.43 
1.55 1.19 
7.51 5.71 

3.58 2.70 

1.66 1.25 

7.59 5.72 
3.43 2.58 
1.53 1.16 
6.75 5.23 
3.00 2.32 
1.35 1.05 



6.84 6.93 
5.35 5.40 

4.06 4.08 
2.96 2.95 

2.07 2.06 
13.95 13.85 

9.15 9.07 
5.84 5.75 
3.64 3.56 
2.22 2.15 

13.31 12.77 
7.84 7.40 
4.53 4.21 

25.52 23.31 

14.07 12.65 
7.64 6.73 
4.07 3.52 

21.42 18.17 
11.10 9.27 
5.74 4.68 

29.08 23.45 
14.80 11.76 

7.51 5.88 
38.14 29.47 
19.31 14.73 

9.91 7.47 



3.35 3.20 

2.64 2.48 

2.00 1.86 
1.47 1.34 
1.04 0.93 
7.13 6.31 
4.72 4.11 

3.01 2.60 
1.87 1.60 
1.13 0.96 
6.67 5.55 
3.81 3.14 
2.13 1.73 

11.58 9.27 

6.12 4.83 

3.13 2.43 
1.56 1.19 
7.55 5.77 
3.60 2.71 
1.66 1.26 
7.59 5.72 
3.42 2.58 
1.53 1.16 
6.77 5.19 
3.01 2.33 
1.35 1.05 



7.88 7.99 
6.13 6.21 
4.65 4.68 
3.36 3.37 

2.34 2.33 
15.61 15.58 
10.11 10.03 

6.38 6.29 
3.92 3.84 

2.35 2.28 
13.85 13.29 

8.05 7.62 

4.59 4.26 
25.58 23.46 
14.07 12.63 

7.59 6.69 

4.03 3.49 
21.22 18.04 
11.06 9.21 

5.68 4.65 
29.05 23.44 
14.80 11.74 

7.51 5.89 
38.13 29.48 
19.38 14.75 

9.90 7.47 



3.90 3.66 

3.05 2.83 
2.32 2.11 
1.69 1.51 
1.19 1.05 

8.10 7.00 
5.29 4.52 

3.34 2.82 

2.06 1.71 
1.22 1.00 
7.04 5.77 
3.98 3.23 
2.19 1.76 

11.71 9.35 

6.12 4.81 

3.11 2.42 
1.54 1.19 
7.51 5.71 

3.57 2.67 
1.65 1.26 

7.58 5.71 
3.43 2.59 
1.53 1.16 
6.78 5.20 
3.02 2.33 

1.35 1.05 
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TABLE XVII. Neutrino flux (m-^sec-^sr-^GeV-^) for -0.6 > cos6'2 > -0.7 



(GeV) 



Kamioka 



Sudbury 



Gran Sasso 



Up Vp 



Norm 



8.31 
6.52 
4.94 
3.57 
2.47 
16.41 
10.52 
6.53 
3.92 
23.09 
13.33 
7.54 
4.20 
22.84 
12.25 
6.43 
3.34 
17.14 
8.74 
4.42 
22.12 
11.08 
55.62 
27.76 
13.89 
70.11 



8.31 4.07 

6.52 3.20 
4.94 2.44 
3.57 1.78 
2.47 1.26 

16.41 8.46 

10.52 5.48 

6.53 3.42 
3.92 2.06 

23.09 12.12 

13.33 6.89 

7.54 3.82 
4.20 2.07 

22.84 10.96 

12.25 5.63 

6.43 2.83 

3.34 1.38 

17.14 6.63 

8.74 3.12 

4.42 1.44 

22.12 6.48 

11.08 2.91 

55.62 12.92 

27.76 5.73 

13.89 2.53 

70.11 11.34 



3.82 
2.97 
2.21 
1.59 
1.10 
7.30 
4.66 
2.88 
1.72 
9.98 
5.60 
3.08 
1.65 
8.64 
4.42 
2.19 
1.06 
5.05 
2.35 
1.09 
4.90 
2.20 
9.90 
4.39 
1.96 
8.86 



5.63 5.71 2.73 2.70 
4.52 4.55 2.19 2.14 
3.52 3.52 1.70 1.63 
2.62 2.61 1.27 1.20 
1.88 1.86 0.92 0.85 

12.96 12.82 6.38 5.82 

8.64 8.50 4.27 3.83 
5.60 5.48 2.77 2.45 
3.54 3.42 1.74 1.51 

21.66 20.84 10.53 9.03 
13.03 12.37 6.22 5.22 

7.66 7.17 3.53 2.95 

4.41 4.06 1.97 1.60 

24.81 22.44 10.59 8.49 

13.67 12.10 5.53 4.35 
7.37 6.40 2.82 2.18 
3.93 3.33 1.38 1.06 

20.56 17.16 6.63 5.03 

10.68 8.72 3.12 2.37 
5.48 4.42 1.44 1.08 

27.95 22.12 6.49 4.90 

14.22 11.09 2.91 2.20 

72.17 55.51 12.93 9.88 

36.55 27.81 5.73 4.39 

18.48 13.91 2.53 1.96 

94.53 70.21 11.32 8.87 



7.40 7.53 3.64 3.49 

5.89 5.96 2.90 2.74 

4.54 4.56 2.23 2.06 

3.33 3.34 1.65 1.50 

2.34 2.34 1.17 1.04 
15.81 15.72 8.02 7.02 
10.31 10.19 5.26 4.53 

6.49 6.38 3.32 2.81 

3.97 3.86 2.02 1.69 

23.68 22.83 11.89 9.81 

13.84 13.17 6.80 5.56 

7.96 7.47 3.77 3.05 

4.51 4.16 2.05 1.64 

25.09 22.67 10.83 8.57 

13.70 12.10 5.59 4.37 

7.35 6.38 2.81 2.18 
3.91 3.31 1.37 1.05 

20.41 17.02 6.59 5.03 

10.67 8.69 3.10 2.35 

5.46 4.41 1.42 1.07 

27.94 22.10 6.47 4.90 

14.20 11.09 2.91 2.21 

72.17 55.50 12.93 9.81 

36.56 27.87 5.74 4.40 

18.51 13.95 2.55 1.96 

94.39 70.16 11.34 8.90 
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TABLE XVIII. Neutrino flux (m-^sec-^sr-^GeV^M for -0.7 > cos 6'., > -0.8 



(GeV) 



Kamioka 



Sudbury 



Vp Up 



Gran Sasso 



Up Up 



Norm 



7.96 
6.32 
4.84 
3.53 
2.45 
16.33 
10.48 
6.48 
3.88 
22.80 
13.08 
7.38 
4.08 
22.17 
11.75 
6.16 
31.88 
16.36 
8.31 
41.84 
21.01 
10.53 
52.76 
26.41 
13.20 
66.52 



7.96 3.89 
6.32 3.11 
4.84 2.37 
3.53 1.75 
2.45 1.23 
16.33 8.29 
10.48 5.36 
6.48 3.32 
3.88 2.00 
22.80 11.62 
13.08 6.54 
7.38 3.61 
4.08 1.94 
22.17 10.17 

11.75 5.18 
6.16 2.57 

31.88 12.42 
16.36 5.88 
8.31 2.74 
41.84 12.58 
21.01 5.64 
10.53 2.51 

52.76 11.18 
26.41 4.92 
13.20 2.18 
66.52 9.75 



3.65 
2.87 
2.16 
1.56 
1.08 
7.17 
4.57 
2.79 
1.66 
9.51 
5.33 
2.91 
1.54 
7.99 
4.03 
1.98 
9.52 
4.46 
2.08 
9.48 
4.27 
1.92 
8.54 
3.82 
1.69 
7.70 



5.33 5.39 
4.33 4.36 
3.41 3.41 
2.57 2.56 
1.86 1.84 
12.92 12.73 
8.65 8.50 
5.62 5.49 
3.53 3.41 
21.61 20.64 
12.88 12.16 
7.55 7.00 
4.32 3.93 
24.30 21.64 
13.32 11.63 
7.16 6.12 
37.98 31.74 
19.90 16.32 
10.32 8.29 
53.00 41.93 
26.96 21.03 
13.70 10.53 
69.47 52.74 
35.19 26.41 
17.73 13.16 
90.59 66.56 



2.57 2.58 
2.08 2.06 

1.63 1.58 
1.23 1.17 
0.89 0.83 
6.23 5.71 
4.19 3.77 
2.71 2.39 
1.69 1.46 

10.14 8.67 
5.94 4.97 
3.36 2.77 

1.85 1.50 
9.82 7.82 
5.05 3.98 
2.53 1.97 

12.33 9.47 

5.86 4.47 
2.74 2.07 

12.40 9.48 

5.64 4.27 
2.51 1.91 

11.16 8.53 
4.92 3.80 
2.18 1.69 
9.79 7.64 



6.59 6.69 3.24 3.11 

5.31 5.36 2.60 2.47 

4.13 4.15 2.01 1.88 
3.07 3.08 1.50 1.37 
2.17 2.17 1.07 0.96 

14.78 14.68 7.34 6.49 

9.71 9.56 4.82 4.19 

6.14 6.03 3.05 2.60 
3.78 3.67 1.86 1.57 

22.70 21.80 10.95 9.12 

13.34 12.62 6.27 5.14 
7.71 7.17 3.49 2.83 
4.39 4.00 1.89 1.51 

24.35 21.75 9.92 7.89 
13.29 11.64 5.11 3.98 

7.14 6.11 2.53 1.96 

37.93 31.68 12.31 9.49 

19.88 16.30 5.83 4.46 

10.29 8.28 2.72 2.06 

53.01 41.86 12.52 9.47 

26.98 21.00 5.64 4.27 

13.70 10.53 2.51 1.92 

69.56 52.79 11.20 8.54 

35.27 26.40 4.92 3.79 

17.77 13.20 2.18 1.69 

90.64 66.58 9.78 7.68 
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TABLE XIX. Neutrino flux (m-^sec-^sr-^GeV-M for -0.8 > cos6l, > -0.9 



(GeV) 



Kamioka 



Sudbury 



Vp Up 



Gran Sasso 



Norm 



7.04 7.04 

5.67 5.67 
4.39 4.39 
3.23 3.23 

22.57 22.57 
15.12 15.12 
9.77 9.77 
6.09 6.09 

3.68 3.68 
21.70 21.70 

12.51 12.51 
7.08 7.08 
3.92 3.92 

21.24 21.24 
11.27 11.27 
5.92 5.92 

30.52 30.52 
15.67 15.67 

7.96 7.96 
40.17 40.17 
20.11 20.11 
10.08 10.08 
50.47 50.47 
25.22 25.22 
12.61 12.61 
63.51 63.51 



3.45 3.25 

2.78 2.57 

2.13 1.95 

1.58 1.42 
11.16 9.84 

7.54 6.57 

4.87 4.19 

3.03 2.59 

1.83 1.53 

10.70 8.82 

6.02 4.94 

3.31 2.69 

1.78 1.43 

9.34 7.34 

4.73 3.70 

2.33 1.80 

11.25 8.57 

5.29 4.02 

2.43 1.86 

11.11 8.41 

4.97 3.77 

2.21 1.69 

9.76 7.51 

4.31 3.33 

1.91 1.50 

8.59 6.78 



6.42 6.46 
5.18 5.20 
4.03 4.02 

2.98 2.96 
20.94 20.78 
14.12 13.96 

9.18 9.02 
5.79 5.65 
3.54 3.42 
21.24 20.30 

12.52 11.79 
7.29 6.72 
4.17 3.76 

23.30 20.58 
12.78 11.05 
6.91 5.83 
36.69 30.32 
19.28 15.60 

9.99 7.93 

51.53 39.95 
26.19 20.13 
13.30 10.08 
67.36 50.47 
34.03 25.29 
17.16 12.62 
87.22 63.48 



3.15 3.01 

2.53 2.38 

1.95 1.80 

1.44 1.31 

10.22 9.13 

6.91 6.07 

4.49 3.89 

2.80 2.40 

1.70 1.43 

9.90 8.28 
5.62 4.65 
3.11 2.55 
1.69 1.36 
8.88 7.09 

4.55 3.59 
2.27 1.77 

11.05 8.49 

5.23 4.00 

2.42 1.84 

11.11 8.36 
4.97 3.77 

2.21 1.68 

9.76 7.49 

4.33 3.33 

1.91 1.50 

8.56 6.75 



6.37 6.43 
5.15 5.19 
4.01 4.02 
2.97 2.96 
20.92 20.79 

14.12 13.99 
9.22 9.10 
5.82 5.69 
3.58 3.46 

21.53 20.54 
12.70 11.95 
7.39 6.83 
4.22 3.82 
23.55 20.75 
12.85 11.12 
6.94 5.84 
36.81 30.32 

19.28 15.59 
10.03 7.91 

51.29 39.84 
26.19 20.11 
13.29 10.09 
67.31 50.42 
34.03 25.26 

17.13 12.60 
87.23 63.52 



3.15 2.97 
2.53 2.36 
1.95 1.79 
1.45 1.30 

10.26 9.11 

6.95 6.08 

4.50 3.90 

2.83 2.42 

1.71 1.44 

10.01 8.35 

5.70 4.72 

3.16 2.58 

1.71 1.38 
9.00 7.16 
4.60 3.62 

2.30 1.78 
11.08 8.50 

5.24 3.99 

2.43 1.84 

11.04 8.36 

4.97 3.77 

2.21 1.68 

9.76 7.51 

4.31 3.36 
1.91 1.49 
8.58 6.73 
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TABLE XX. Neutrino flux (m-^sec-^sr-^GeV-^) for -0.9 > cos6'^ > -1.0 



(GeV) 



Kamioka 



Sudbury 



Gran Sasso 



Up Vp 



Norm 



6.78 6.78 
5.51 5.51 
4.30 4.30 
3.18 3.18 
22.38 22.38 
15.05 15.05 
9.74 9.74 
6.07 6.07 
3.65 3.65 
21.51 21.51 
12.36 12.36 
6.95 6.95 
3.84 3.84 
20.67 20.67 
10.95 10.95 
5.71 5.71 
29.41 29.41 
15.04 15.04 
7.63 7.63 
38.59 38.59 
19.33 19.33 
9.70 9.70 
48.53 48.53 
24.28 24.28 
12.08 12.08 
60.78 60.78 



3.30 3.15 

2.66 2.51 
2.06 1.91 
1.53 1.39 

10.84 9.68 
7.34 6.44 
4.75 4.11 
2.95 2.52 
1.77 1.49 

10.22 8.50 
5.77 4.73 
3.15 2.55 
1.68 1.34 

8.67 6.86 
4.40 3.41 
2.14 1.65 

10.21 7.80 
4.77 3.63 
2.19 1.67 
9.91 7.52 
4.43 3.36 
1.97 1.50 

8.66 6.70 
3.84 2.99 
1.70 1.34 

7.67 6.06 



8.39 8.47 

6.76 6.80 
5.22 5.22 
3.81 3.79 

26.28 26.09 
17.32 17.13 
10.98 10.82 
6.73 6.57 
4.02 3.87 

23.48 22.26 

13.49 12.59 
7.67 7.01 
4.29 3.84 

23.70 20.66 
12.80 10.91 
6.85 5.68 
36.14 29.39 
18.84 15.03 

9.77 7.65 
50.21 38.57 

25.51 19.36 
12.94 9.70 

65.52 48.58 
33.06 24.31 
16.61 12.10 
84.31 60.81 



4.17 3.84 
3.34 3.03 
2.55 2.27 
1.86 1.63 

12.91 11.13 

8.55 7.25 

5.40 4.52 

3.28 2.71 

1.91 1.57 

10.84 8.82 

5.97 4.82 

3.21 2.57 

1.69 1.34 

8.68 6.83 
4.36 3.39 
2.13 1.64 

10.18 7.79 

4.75 3.62 

2.18 1.66 

9.84 7.56 
4.43 3.37 
1.97 1.50 

8.70 6.71 

3.85 2.99 

1.69 1.34 
7.66 6.08 



7.42 7.49 
6.02 6.04 
4.67 4.67 
3.44 3.42 

23.92 23.74 

15.93 15.75 
10.24 10.06 

6.36 6.20 
3.84 3.69 
22.63 21.52 
13.16 12.29 
7.54 6.89 
4.26 3.79 
23.55 20.51 
12.80 10.90 
6.84 5.68 
36.20 29.35 
18.85 15.03 
9.79 7.65 
50.19 38.73 
25.50 19.35 
12.92 9.71 
65.41 48.57 
32.98 24.26 
16.60 12.10 
84.26 60.81 



3.67 3.44 

2.96 2.72 

2.27 2.05 

1.67 1.48 

11.68 10.19 

7.78 6.72 

4.96 4.23 
3.04 2.56 
1.81 1.50 

10.31 8.52 

5.75 4.70 

3.13 2.53 

1.66 1.33 

8.60 6.79 

4.33 3.39 

2.13 1.64 

10.17 7.79 

4.75 3.64 

2.20 1.66 

9.90 7.52 

4.43 3.37 

1.97 1.50 
8.70 6.71 
3.84 2.97 
1.70 1.33 
7.66 6.06 
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